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Preface 


This  investigation  brings  the  APTT  Low  Speed  Shock  Tube  back  into  use  on  a  study  of  the 
effects  of  transpiration  cooling  on  the  boundary  layer  in  a  Mach  2.0  nozzle.  The  experiment  was 
performed  in  parallel  with  a  heat  transfer  reduction  study  carried  out  using  the  same  equipment. 
The  focus  of  the  study  was  characterizing  the  effects  that  the  increased  boundary  layer  thickness 
had  on  the  performance  of  the  nozzle  in  terms  of  boundary  layer  displacement  thickness  and  exit 
Mach  number.  This  is  the  first  in  a  new  series  of  experiments  that  are  to  be  run  on  this  nozzle, 
and  a  great  amount  of  headway  has  been  made  during  the  course  of  this  project.  That  headway 
would  not  have  been  possible,  however,  without  the  help  of  several  people. 

Many  thanks  go  out  to  my  advisor.  Dr.  Rodney  Bowersox,  and  Lt  Col  Jerry  Bowman  for 
their  knowledge  of  the  subject,  the  wide  berth  allowed  me  to  perform  my  work,  and  answers  to  aU 
my  questions.  Thanks  and  appreciation  go  to  Mr  .Tim  Hancock  of  the  AFTT  Machine  Shop  for  his 
excellent  work  on  the  model.  He  was  asked  to  do  many  unfamiliar  things  and  did  so  with 
professionalism,  creativity,  and  precision.  Also,  Mr.  Andy  Pitts  and  the  other  members  of  the 
AFTT/ENY  Laboratory  Staff  receive  thanks  for  helping  with  all  the  hardware  concerns 
surrounding  this  project.  Special  thanks  goes  to  Andy  for  his  help  in  setting  up  the  timing  circuit 
that  allowed  me  great  success  at  taking  shadowgraphs. 

I  enjoyed  my  time  working  in  Rm.l46,  and  will  take  all  the  things  I  have  learned  from  all 
of  you  and  hopefully  put  them  to  good  use  at  my  next  assignment  and  beyond. 
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Abstract 

The  effects  of  transpiration  cooling  through  a  porous  wall  on  nozzle  performance  was 
investigated.  The  experiments  were  performed  in  the  AFIT  Low  Pressure  shock  tunnel.  The 
tunnel  was  fitted  with  a  two  dimensional,  Mach  2.0  characteristic  nozzle  with  one  wall 
constructed  of  porous  material.  Blowing  ratios  from  0.4%-1.2%  of  the  mainstream  flow  were 
studied.  This  study  extends  current  flat  plate  knowledge  of  transpiration  cooling  by  including 
compressibility  and  a  severe  favorable  pressure  gradient. 

Boundary  layer  growth  as  a  result  of  the  transpiration  was  measured  using  flow 
visualization,  and  the  effect  of  that  growth  on  the  exit  Mach  number  profile  was  investigated. 
Analysis  of  experimental  data  showed  increases  of  up  to  47%  in  boundary  layer  thickness  with 
transpiration,  and  because  the  blowing  ratios  used  in  this  experiment  were  very  low  due  to 
equipment  constraints,  the  resulting  trends  show  no  sign  of  leveling  off  leaving  much  room  for 
continued  research  of  this  kind  at  higher  blowing  ratios.  The  results  do,  however,  show  that 
performance  losses  due  to  transpiration  are  such  that  transpiration  may  be  a  viable  nozzle  cooling 
method  even  up  to  high  blowing  ratios  without  having  large  adverse  effects  on  performance. 

This  study,  taken  together,  with  a  parallel  study  on  heat  transfer  reduction  using  the  same 
equipment  endeavors  to  begin  to  characterize  the  tradeoff  between  the  performance  gain  due  to 
heat  transfer  reduction  and  the  loss  due  to  boundary  layer  growth  in  a  transpiration  cooled  nozzle. 
Defining  the  optimal  blowing  ratio  to  minimize  heat  transfer  and  optimize  performance  would  be 
the  final  goal  of  such  an  effort. 
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INVESTIGATION  OF  BOUNDARY  LAYER  AND 
PERFORMANCE  EFFECTS  OF  TRANSPIRATION  COOLING 
THROUGH  A  POROUS  PLATE  ON  A  ROCKET  NOZZLE 

I.  Introduction 


1.1  Background 

Currently  there  are  quality  rocket  engine  systems  and  high  temperature  wind  tunnel 
facihties  in  use.  Further  improvements  in  the  performance  of  these  systems  requires  higher 
operating  temperatures.  A  problem  is  that,  without  cooling,  operating  at  higher  temperatures 
tends  to  compromise  the  structural  integrity  of  the  thrust  chamber.  For  example,  in  hypersonic 
wind  tunnels,  for  example,  erosion  of  the  tunnel  wall  sometimes  occurs  due  to  high  operating 
temperatures.  Erosion  can  cause  debris  to  travel  down  the  tunnel  that  often  breaks  expensive  test 
articles  because  of  impact  speeds.  Nozzle  heat  transfer  and  material  thermal  limitations,  then, 
remain  a  hmiting  factor  in  the  performance  of  modem  rocket  engines.  Therefore,  heat  transfer 
reduction  is  the  key  to  improving  the  efficiency  stmctural  integrity  of  rocket  systems,  but  the 
relationship  between  heat  transfer  and  the  boundary  layer  present  in  the  flow  field  is  also  very 
important. 

Heat  transfer  and  boundary  layer  thickness  in  a  nozzle  are  inversely  related  phenomena 
(Hill  and  Peterson,  1992:544).  The  boundary  layer  thickness  has  been  demonstrated  analytically 
by  Bartz  (1955:1243)  to  be  a  minimum  at  the  throat  of  a  supersonic  nozzle.  In  addition,  its 
thickness  at  the  throat  is  not  highly  dependent  on  its  thickness  before  the  throat  or  the  rate  of 
convergence  of  the  nozzle  after  the  throat  (Hill  and  Peterson,  1992:545).  This  means  that  in 
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virtually  every  case,  the  boundary  layer  thickness  will  be  a  minimum  at  the  throat  of  a  nozzle,  and 
this  is  a  very  important  point  in  relation  to  heat  transfer. 

A  turbulent  boundary  layer  provides  insulation  to  nozzle  walls  that  makes  heat  transfer  to 
the  walls  less  than  if  no  boundary  layer  were  present.  The  thicker  the  boundary  layer,  the  more 
insulation  it  provides,  and  the  lower  the  heat  transfer  through  it.  This  is  the  inverse  relationship 
mentioned  above,  and  it  implies  that  heat  transfer  will  be  greatest  in  a  rocket  nozzle  at  the  throat 
because  the  boundary  layer  is  a  minimum  there  (Hill  and  Peterson,  1992:545).  A  further 
implication  is  that  heat  transfer  will  be  a  maximum  at  the  throat  in  nearly  every  case  because  of 
throat  boundary  layer  thickness  being  minimum  in  nearly  every  case.  This  relationship  has  been 
validated  both  analytically  and  experimentally  (Hill  and  Peterson,1992:551).  Actively  cooled 
rocket  nozzles  allow  the  high  combustion  temperatures  necessary  for  high  performance  while 
maintaining  the  structural  integrity  of  the  nozzle.  Several  methods  of  active  cooling  have  been 
employed  to  address  this  problem,  including  regenerative,  film,  and  transpiration  cooling. 

Regenerative  cooling  involves  pumping  a  liquid  through  channels  surrounding  the  outside 
of  the  combustion  chamber.  This  method  has  been  used  extensively  because  of  its  relative 
simplicity  and  low  cost.  Also,  the  heat  absorbed  by  the  coolant  is  not  wasted  because  it  augments 
the  initial  energy  content  of  the  propellant  prior  to  injection,  increasing  the  exhaust  velocity 
slightly  (Sutton,  1992:290).  One  drawback  to  regenerative  cooling  is  the  additional  pumping 
requires  to  force  coolant  through  the  coolant  lines.  This  reduces  the  turbopump’s  ability  to 
pressurize  the  combustion  chamber  and  leads  to  lower  expansion  ratios  and  associated  lower 
specific  impulses. 
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Film  cooling  involves  injecting  low  temperature  gas  or  liquid  through  one  or  several 
discrete  holes  in  the  nozzle  wall  to  establish  a  protective  film  on  the  surface.  Because  the  coolant 
is  inside  the  thrust  chamber,  it  accepts  heat  directly,  which  makes  it  a  more  effective  cooling 
method  than  regenerative  cooling.  It  is  relatively  easy  to  implement  and  has  been  used  in  thrust 
chambers  for  many  years  (Sutton,  1992:294).  The  drawback  of  film  cooling  is  that  it  requires 
large  injection  mass  flow  rate  per  unit  area  (10-300%  mainstream  flow),  which  can  cause  thrust 
losses  due  to  disturbance  of  the  primary  flow  (Azevedo,  1993:44). 

Transpiration  cooling  is  essentially  the  limiting  case  of  film  cooling  because  it  involves 
pushing  gas  or  liquid  uniformly  through  an  area  of  porous  wall  material  (Sutton,  1992:294).  It  can 
be  thought  of  as  an  infinite  number  of  film  cooling  ports  with  zero  distance  between  them. 
Though  the  concept  is  similar  to  film  cooling,  the  effects  on  heat  transfer  and  boundary  layer  are 
very  different  because  the  fluid  enters  the  flow  through  the  porous  material  at  very  low  velocity 
and  does  not  cause  as  great  a  disturbance  to  the  flow  field.  This  cooling  technique  has  been 
successfully  used  for  cooling  injector  faces  on  the  moon  launch  vehicle  and  the  Space  Shuttle 
Main  Engines,  but  it  has  not  been  used  in  cooling  the  thrust  chamber  or  nozzle  regions  of  large 
rocket  engines.  Due  to  the  steep  pressure  gradients  along  the  inner  wall  of  the  nozzle,  especially 
near  the  throat,  proper  cooling  requires  a  variable  porosity  and/or  thickness  wall  material.  The 
manufacture  of  large,  complex  shapes  of  porous  materials  is  a  difficult  problem,  and  has  been  a 
great  challenge  to  transpiration  cooling  in  this  region  (Sutton,  1992:294). 

As  mentioned  above,  transpiration  cooling  has  been  used  very  little  in  practice,  and  no 
generalized  analytical  relations  exist  for  the  prediction  of  heat  flux  or  boundary  layer  growth  in 
configurations  with  pressure  gradient,  as  in  the  case  of  supersonic  nozzles  (Beitel,  1993:49).  It  is 
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possible  that  heat  transfer  rate  reduction  similar  to  or  better  than  film  cooling  can  be  realized  with 
decreased  flow  disturbances,  indicating  that  transpiration  cooling  could  be  a  more  attractive 
method  than  film  or  regenerative  cooling.  Therefore,  this  investigation  proposes  to  study  the 
effects  of  transpiration  cooling  in  a  supersonic  nozzle  on  boundary  layer  growth  and  performance 
at  low  blowing  ratios. 

1.2  Problem 

Transpiration  cooling  has  a  significant  effect  on  the  turbulent  boundary  layer  structure 
while  being  a  very  effective  way  to  protect  surface  from  hot  mainstream  flow  (Kays  and 
Crawford,  1980: 179,223).  Semiempirical  relationships  have  been  developed  for  predicting  the 
effect  of  transpiration  on  turbulent  boundary  layers  over  a  flat  plate,  and  additional  data  has  been 
taken  to  validate  these  relations.  As  mentioned  above,  transpiration  cooling  has  been  used  very 
little  in  practice,  and  only  practical  data  will  show  its  true  benefits.  Heat  transfer  rate  reduction 
with  transpiration  cooling  may  be  similar  to  film  cooling,  and  because  of  decreased  flow 
disturbances,  reduced  performance  loss  for  the  same  heat  transfer  reduction  could  be  achieved. 
Therefore,  transpiration  cooling  could  be  an  attractive  alternative  to  film,  but  practical  data  is 
needed  to  show  this  benefit. 

The  problem  addressed  here,  then,  was  the  lack  of  knowledge  of  the  effects  of 
transpiration  cooling  in  supersonic  nozzles  because  data  of  this  type  was  needed  to  show  how 
transpiration  cooling  could  improve  current  rocket  systems.  Specific  parameters  that  were 
investigated  include  boundary  layer  thickness,  exit  Mach  number,  and  specific  impulse,  all  of 
which  affect  the  performance  of  the  nozzle. 
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1.3  Summary  of  Current  Knowledge 

Since  the  goal  of  this  research  was  showing  the  effect  of  transpiration  on  boundary  layer 
thickness  and  performance,  previously  collected  knowledge  in  each  category  was  desired  for 
comparison.  The  situation  here  involved  a  severe  favorable  pressure  gradient  (i.e.,  a  supersonic 
nozzle),  and  it  appeared  that,  as  Azevedo  asserts  (1993:43),  most  studies  up  to  this  time  had 
concentrated  on  the  heat  transfer  aspect  of  transpiration  cooling.  In  addition,  most  of  those 
studies  used  flat  plates.  There  was  not  a  wealth  of  data  available  concerning  the  effects  of 
transpiration  cooling  on  boundary  layer  growth  and  even  less  concerning  its  effect  on  nozzle 
performance,  but  following  is  a  synopsis  of  the  information  collected. 

While  transpiration  cooling  has  never  made  it  into  wide  practice,  it  is  not  a  new  idea. 
Since  the  1950s,  it  has  been  the  subject  of  research  aimed  at  cooling  aerodynamic  surfaces  subject 
to  high  heating  (Kays  and  Moffat,  1975:224).  Aerojet  has  dealt  in  systems  employing  porous 
surfaces  since  1946,  and  in  1967  they  fired  the  Aerojet  ARES  to  demonstrate  the  technology. 
This  system  produced  445  kN  (100,000  Ibf)  of  thrust  with  a  20.7  MPa  (3000  psi)  transpiration 
cooled  thrust  chamber.  Aside  from  this,  however,  transpiration  cooled  thrust  chambers  are 
virtually  foreign  to  the  commercial  rocket  industry.  Research  in  the  area  has  continued,  and  a 
group  at  Stanford  University  collected  a  large  amount  of  data  between  1965-74  that  makes  up 
much  of  the  transpiration  cooling  knowledge  base  (Kays  and  Moffat,1975:223). 

For  instance,  Simpson  accomplished  a  study  of  the  boundary  layer  on  a  porous  wall  with 
suction  and  blowing  in  1968.  His  study  showed  boundary  layer  growth  in  both  cases 
(Schetz,  1984:205),  and  plots  from  his  results  are  shown  in  Figure  1.1.  In  contrast,  Antonia  and 
Fulachier  noticed  a  23%  decrease  in  boundary  layer  thickness  in  the  presence  of  suction 
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(1989:431).  Scott  and  others  made  similar  measurements  by  blowing  helium  through  a  porous 
plate  and  also  saw  boundary  layer  growth  as  a  result  (Schetz,1984:273),  although  it  was  not  as 
pronounced  as  in  the  Simpson  case  (See  Figure  1.2).  In  Figures  1.1  and  1.2,  5*  represents 
boundary  layer  displacement  thickness,  and  actual  data  points  are  represented  by  the  filled 
squares. 


(a)  -  Suction  (b)  -  Blowing 

Figure  1.1  -  Boundary  Layer  Growth  w/  Blowing  (Schetz, 1984:205) 


Figure  1.2  -  Boundary  Layer  Thickness  w/  Helium  Injection  (Schetz, 1984:276) 
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.  d  fcm) 


The  lines  accompanying  the  data  in  Figures  1.1  and  1.2  represent  the  results  of  numerical 
methods  developed  by  Schetz  and  Favin  for  predicting  the  boundary  layer  thickness  with  injection 
or  suction  (Schetz, 1984:202,273).  Their  analysis  used  the  Reichardt  eddy  viscosity  model 
extended  to  injection  and  applied  it  to  the  above  cases  for  their  predictions  (Schetz,  1984:202). 
Schetz  and  Favin’s  method  predicted  very  accurately  the  experimental  data  collected  in  Simpson’s 
suction  case  (Schetz,  1984:205)  and  Scott’s  blowing  case  (Schetz,  1984:276),  but  its  predictions 
were  slightly  low  in  the  Simpson  blowing  case  (Schetz,  1984:205). 

The  contrasting  results  of  Simpson  versus  Antonia  and  Fulachier  and  the  varying  results  of 
Schetz  and  Favin’s  prediction  methods  point  out  the  complexity  of  a  transpired  turbulent 
boundary  layer.  Another  result  from  Schetz  (White,  199 1:435)  states  that  a  porous  wall  actually 
increases  skin  friction  (related  to  Cf)  along  a  nozzle  wall  which  would  decrease  the  boundary  layer 
thickness.  This  also  means  that  a  certain  amount  of  blowing  would  be  required  to  counteract  the 
increase  in  Cf  due  to  the  porous  plate  before  any  cooling  can  be  achieved  (White,  199 1:436).  In 
addition  to  these  effects,  roughness  of  a  porous  surface  can  also  affect  boundary  layer  growth. 
Therefore,  the  conflicting  results  above  are  not  disturbing  because  the  interrelationships  between 
suction,  blowing,  and  roughness  are  not  well  understood.  Figures  1.1  and  1.2  only  show  that  it  is 
very  difficult  to  predict  the  exact  effect  transpiration  cooling  will  have  on  a  given  boundary  layer. 

An  older  study  by  Goldstein,  Shavit,  and  Chen,  involved  injecting  air  into  air  through  a 
section  of  porous  material  on  a  flat  plate.  They  measured  a  rapid  increase  in  boundary  layer 
thickness  at  the  beginning  of  the  porous  material  with  slower  increases  further  downstream 
(Figure  1.3).  Also,  the  thickness  reached  a  definite  maximum  at  a  point  downstream  of  the 
injection  point,  and  far  downstream  the  boundary  layer  resumed  the  normal  shape  of  a  boundary 
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layer  on  a  flat  plate  (Goldstein,  Shavit,  and  Chen,1965:355).  The  results  of  the  Goldstein  et.  al. 
study  provided  an  indication  of  what  to  expect  in  this  experiment  as  a  similar  situation  existed. 
The  introduction  of  a  pressure  gradient  in  this  researeh  did  change  the  situation,  and  the  results 
seen  here  were  quite  different. 


Figure  1.3  -  Flat  Plate  Boundary  Layer  Growth  w/  Injection  (Goldstein  et.al., 1965:356) 

The  Goldstein  et.  al.  study  also  points  out  that  when  comparing  injection  ratios  it  may  be 
better  to  compare  injection  ratio  per  unit  length  than  just  plain  injection  ratio  (Goldstein  et 
al,  1965:355).  This  may  provide  a  method  of  relating  the  results  of  experiments  using  different 
geometries  to  each  other. 

One  performance  study,  accomplished  by  Azevedo,  used  three  rows  of  film  cooling  holes 
in  a  diverging  nozzle  for  injection  and  measured  the  effect  of  blowing  on  the  thrust  efficiency  of 
the  nozzle.  The  parameters  investigated  in  this  study  included  injection  angle  (15°-100°),  nozzle 
divergence  angle  (1°-10°),  and  blowing  ratio  (0.01-0.08).  Because  transpiration  cooling  is 
essentially  the  limiting  case  of  90°  film  cooling,  only  data  for  the  100°  injection  angle  case  was  of 
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interest  because  it  was  the  closest  to  the  transpiration  condition.  Thrust  efficiency  (Cv)  was 
referenced  to  the  ideal  thrust  the  nozzle  would  produce  under  isentropic  conditions 
(Azevedo,  1993:47),  and  results  were  reported  in  terms  of  ACv  from  the  non  blowing  to  the 
blowing  case. 

Not  surprisingly,  Azevedo  found  that  injection  at  100°  caused  the  greatest  ACyi  more  than 
-0.06  or  a  performance  loss  of  6%  in  terms  of  thrust  (Azevedo,  1993:48).  Some  of  Azevedo’s 
results  are  shown  in  Figures  1.4a  and  b.  The  largest  injection  angle  caused  the  largest 
performance  loss  because  larger  injection  angles  cause  greater  flow  disturbances,  which  take  more 
energy  to  overcome. 


A 


(a)  (b) 

Figure  1.4  -  Thrust  Losses  w/  Film  Cooling  in  a  Nozzle  (Azevedo, 1993:48) 


Azevedo’s  study  applied  to  this  experiment  because  it  gave  an  idea  of  the  effect  that 
injection  by  transpiration  could  have  on  the  performance  of  a  nozzle.  However,  Azevedo’s  results 
were  very  specific  to  the  geometries  he  tested  because  it  is  very  difficult  to  translate  the  results  of 
one  experiment  involving  a  pressure  gradient  to  another.  Therefore,  Azevedo’s  results  probably 
had  little  bearing  on  this  experiment  except  as  an  order  of  magnitude  comparison.  It  was 
expected  that  results  in  this  experiment  would  show  smaller  performance  losses  than  were 
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reported  by  Azevedo  because  the  injection  velocities  were  much  smaller,  which  caused  smaller 
flow  disturbances. 

In  another  performance  study  by  Aerojet,  a  platelet  design  for  a  porous  wall  was  used. 
This  type  of  design  addressed  the  problem  of  varying  blowing  ratios  due  to  varying  pressure  drops 
along  the  porous  wall  by  decreasing  port  size  with  increasing  pressure  drop.  Flow  through  a 
porous  wall  of  uniform  porosity  wall  always  seeks  the  path  of  least  resistance,  which  is  where  the 
greatest  pressure  drop  occurs.  In  a  nozzle,  the  pressure  drop  increases  along  the  nozzle  while 
temperature  decreases.  Therefore,  flow  through  a  nozzle  wall  of  uniformly  porous  material  is  not 
directed  to  the  throat,  which  is  the  hot  spot.  A  platelet  design  better  directs  the  cooling  flow  to 
hot  spots  in  a  nozzle  instead  of  allowing  flow  to  go  to  low  temperature  areas,  which  do  not 
require  as  much  cooling. 

The  Aerojet  study  was  a  combustion  study  that  used  transpiration  cooling  with  RP-1 
liquid  rocket  fuel  as  the  coolant.  It  employed  injection  flow  rates  from  9.98-61.2  kg/s  (22.0- 
135.0  Ibm/s)  at  chamber  pressures  from  8.2-19.3  MPa  (1192-2804  psi)  and  compared  the  engine’s 
performance  to  an  engine  using  regenerative  cooling.  Besides  reducing  the  heat  transfer  by  an 
impressive  97%,  with  transpiration  cooling  they  were  able  to  actually  increase  the  rocket’s 
specific  impulse  by  as  much  as  10  sec  over  regeneratively  cooled  systems  (May  and 
Burkhardt,  199 1:53).  Since  specific  impulse,  or  Isp,  is  a  parameter  of  interest  in  this  experiment, 
this  was  an  interesting  result. 
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1.4  Objectives  and  Scope 


In  one  sense,  using  transpiration  cooling  in  a  model  with  complex  geometry  and  a  pressure 
gradient  made  the  scope  of  this  experiment  broad  because  it  pushed  the  envelope  of  current 
knowledge.  At  the  same  time,  however,  the  scope  was  necessarily  small  for  the  same  reason. 
Extending  the  envelope  of  current  knowledge  opened  up  nearly  endless  permutations  of 
combinations  of  porous  material  thickness,  porosity,  area,  position  in  the  nozzle,  and  blowing 
ratio.  This  research  only  hoped  to  scratch  the  surface  of  new  possibilities  with  transpiration 
cooling  and  begin  relating  the  parameters  of  interest  to  each  other.  Test  conditions  were  kept  as 
constant  as  possible  over  the  many  test  runs  conducted,  and  data  was  collected  using  a  single 
nozzle  geometry  with  a  single  thickness  (made  up  of  four  layers  of  equal  thickness  material)  and 
area  of  porous  material  over  a  range  of  low  blowing  ratios  (0.0  -  0.01 1). 

Within  the  above  scope,  the  objectives  of  this  research  were  to  understand  the  differences 
between  flow  over  a  porous  and  non  porous  wall  with  no  blowing,  to  relate  blowing  ratio  to 
boundary  layer  thickness  on  a  porous  wall,  to  relate  blowing  ratio  to  the  uniformity  and 
magnitude  of  the  exit  Mach  number  across  the  exit  plane  of  the  nozzle,  and  relate  the  above 
results  to  the  performance  of  the  nozzle. 

1 .5  Methodology 

A  Mach  2.0  characteristic  nozzle  was  designed  and  installed  in  the  AFTT  low  pressure 
shock  tunnel,  which  can  produce  high  pressure,  high  temperature  chamber  conditions.  Chamber 
pressure  and  temperature  were  approximately  482  kPa  (70  psi)  and  475  K  respectively.  One  side 
of  the  nozzle  (the  non  blowing  side)  was  constructed  with  no  modifications  as  a  control  in  the 
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experiment,  and  the  other  (the  blowing  side)  was  modified  to  accommodate  transpiration  by 
replacing  much  of  the  wall  with  porous  material. 

The  run  times  in  the  shock  tunnel  were  very  short  (i.e.,  milliseconds)  so  high  speed, 
precision  transducers  were  required  for  data  collection.  A  total  of  1 1  pressure  transducers  were 
used  to  collect  the  data  necessary  for  complete  analysis  of  transpiration  cooling  effects. 
Transducers  were  placed  along  both  walls  of  the  supersonic  portion  of  the  nozzle,  in  the 
converging  section  of  the  nozzle,  in  the  constant  pressure  blowing  plenum,  in  a  pitot  probe  for 
exit  pressure  measurements,  and  along  the  top  of  the  tunnel  for  shock  speed  measurements. 

Five  data  sets  were  collected  during  this  experiment:  one  baseline  set  and  four  with 
blowing.  Many  runs  with  no  blowing  in  the  nozzle  were  conducted  to  refine  testing  procedures, 
test  conditions,  and  completely  understand  the  model’s  flow  features  at  baseline  conditions  in 
terms  of  chamber  conditions,  wall  Mach  number,  and  exit  Mach  number  profile.  Shock  speed 
was  also  determined  to  have  a  large  effect  on  chamber  conditions  during  this  phase  of  the 
experiment.  Chamber  conditions,  wall  Mach  number,  and  exit  Mach  number  data  were  also 
taken  with  blowing  for  comparison  with  baseline  data.  As  mentioned  earlier,  only  low  blowing 
ratios  were  considered. 

Another  source  of  data  besides  the  transducers  was  shadowgraph  photography 
accomplished  during  this  experiment  because  they  were  the  only  source.  Shadowgraphs  were 
used  during  the  baseline  phase  of  the  experiment  to  visualize  the  initial  shock  wave  propagation 
through  the  nozzle,  initial  start  up  dynamics,  established  flow,  and  unstart  dynamics.  These 
photographs  were  also  used  to  estimate  the  run  time  of  the  nozzle.  In  addition,  the  effects  that 
the  porous  wall  without  blowing  had  on  the  flow  field  were  determined  from  the  photographs. 
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Boundary  layer  thickness  measurements  were  taken  from  shadowgraphs  of  established  flow  at 
each  blowing  ratio  to  see  the  effect  of  transpiration  on  the  boundary  layer  thickness.  Finally, 
shadowgraphs  proved  useful  at  other  times  for  solving  data  collection  problems. 

Transducer  data  for  each  run  was  reduced  using  computer  programs  written  specifically 
for  this  experiment.  Run  times  varied  slightly  for  each  run,  which  did  cause  variations  in  the 
results.  The  calculated  average  values  for  data  presented  here  include  a  Student  t  99% 
confidence  interval. 

Blowing  ratio  was  defined  throughout  this  experiment  as: 


(1-1) 


where;  B  =  blowing  ratio 

s  =  aic  distance  along  blowing  area 
p  =  gas  density 
u  =  gas  velocity 
i  =  injected  gas  condition 
oc  =  freestream  gas  condition 


This  definition  for  blowing  ratio  was  chosen  because  it  removed  area  dependence  from 
the  calculation,  and  it  took  into  account  the  recommendation  of  Goldstein  et.  al.  (1965:355)  that 
using  a  blowing  ratio  per  unit  length  is  better  than  just  blowing  ratio  alone. 
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II.  Theory 


Before  embarking  on  the  study  of  transpiration  effects  on  boundary  layer  thickness  and 
performance,  some  boundary  layer  theories  related  to  this  study  should  be  reviewed.  Specifically, 
this  study  was  interested  in  the  effects  of  boundary  layer  growth  in  the  presence  of  a  favorable 
pressure  gradient  on  the  performance  of  a  nozzle.  Therefore,  only  theories  leading  to  boundary 
layer  thickness  calculations  are  reviewed  here.  First,  however,  some  facts  and  basic  boundary 
layer  relations  should  be  reviewed. 

2.1  Turbulent  Boundary  Layer 

Prandtl  first  introduced  the  concept  of  a  viscous  boundary  layer  in  1904  and  showed  that  a 
flow  field  could  be  broken  into  a  thin  viscous  layer  near  the  wall  and  the  freestream,  which  is 
essentially  nonviscous  (Hill  and  Peterson,1992:95).  Since  then,  knowledge  of  the  boundary  layer 
has  flourished,  and  distinctions  between  laminar  and  turbulent  layers  have  been  made  as  well  as 
characterization  of  the  transition  from  one  type  to  the  other.  This  review  will  concentrate  on 
turbulent  boundary  layer  theory  because  this  type  of  boundary  layer  was  expected  in  this 
experiment. 

Engineering  applications  such  as  airfoils,  inlets,  and  turbine  blades  that  can  have  a 
turbulent  boundary  layer  almost  always  perform  better  because  of  the  superior  pressure  recovery 
that  can  be  obtained  with  turbulent  over  laminar  boundary  layers.  Therefore,  it  is  generally 
advantageous  to  have  a  turbulent  boundary  layer.  Turbulence  is  a  very  complex  phenomenon, 
and  completely  analyzing  it  is  a  nearly  impossible  task  (White,  199 1:394).  There  is  currently  no 
purely  analytical  treatment  of  turbulence,  and  therefore,  understanding  it  depends  on  experimental 
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observation  (Hill  and  Peterson,  1992: 111).  Equations  have  been  developed  semiempirically  to 
deal  with  turbulent  boundary  layer  analysis,  but  as  more  data  was  collected,  the  basic  boundary 
layer  equations  were  extended.  These  extensions  provided  tools  for  dealing  with  turbulent 
boundary  layers  in  a  variety  of  situations  including  the  transpired  boundary,  which  was  of  interest 
here. 

The  equations  given  in  the  following  development  assume  a  steady,  two-dimensional, 
incompressible  boundary  layer.  Also,  this  development  can  be  found  in  a  more  complete  form  in 
Hill  and  Peterson  (1992:102-105).  The  development  begins,  as  always,  with  a  control  volume 
(Fig  2.1)  and  the  continuity  and  momentum  equations,  Eq.2-1.  Velocities  in  the  x  and  y  direction 
are  represented  by  u  and  v  respectively. 


Figure  2.1  -  Control  Volume 
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Integrating  the  momentum  across  the  boundary  layer  is  a  useful  approximation  for  solving 
the  boundary  layer  equations,  and  Eq.2-2  shows  the  general  form  of  the  momentum  integral 
equation  after  algebraic  manipulation  from  its  original  form  in  Eq.2-lb 

\^[u{U -u)]dy  +  ^\{U -u)dy=+—  (2-2) 

0  o  X  ax  Q  p 

where:  =  an  indefinitely  large  quantity 

To  =  wall  shear  stress 

p  =  flow  density 

u  =  velocity  of  flow  inside  the  boundary  layer 
U  =  freestream  velocity 

The  quantity  (U  -  u)  is  called  the  mass  flow  defect,  which  is  a  result  of  the  decreased  flow 
velocity  in  the  boundary  layer  due  to  friction.  In  this  region  of  the  boundary  layer,  the  mass  flow 
rate  is  actually  less  than  the  freestream  flow.  For  a  nozzle  to  enclose  only  uniform  velocity 
freestream  flow,  the  boundary  layer  would  have  to  remain  in  place,  and  the  walls  bounding  the 
flow  would  have  to  be  displaced  inward.  The  flow  region  excluded  by  moving  the  walls  inward  is 
called  the  displacement  thickness  (5*)  because  this  region  is  unavailable  to  the  freestream. 
Displacement  thickness  is  defined  by  Equation  2-3a.  The  quantity  u(U  -  u)  is  similarly  defined  as 
the  momentum  defect,  and  has  an  associated  thickness  (0)  defined  by  Eq.2-3b.  The  upper  limits 
of  integration  can  be  infinity  in  Eq.  2-3  because  outside  the  boundary  layer  u=  U,  and  the  integral 
goes  to  zero. 

(2-3a) 

0  ^ 


16 


(2-3b) 


where:  5*  =  displacement  thickness 
6  =  momentum  thickness 


«  1  M  M  . 

0  =  1  — (1— — )dy 
^  U  U 


Plugging  these  definitions  into  Eq.2-2  gives  the  most  commonly  used  forms  of  the 
momentum  integral  equation,  Eq.2-4.  This  equation  applies  to  both  laminar  and  turbulent 
boundary  layers,  but  requires  knowledge  of  the  relationship  between  5*  and  6  to  solve.  This 
knowledge  can  be  gained  with  use  of  a  shape  function  (not  considered  here)  or  by  numerical 
momentum  integral  methods. 


dU 

dx 


(ud*)- 


(2-4) 


or 


dQ  e  dU 

— +(2  +  //)— — 


where:  H=Q  /d* 

U  =  freestream  velocity 
Cf=  coefficient  of  friction 


The  general  boundary  layer  equations  developed  above  have  been  extended  and  applied  to 
turbulent  boundary  layers  with  a  pressure  gradient.  Solutions  to  the  equations  for  this  case  can 
also  be  solved  by  at  least  two  numerical  momentum  integral  methods,  which  are  discussed  here. 
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First,  a  formulation  by  Moses  estimates  the  boundary  layer  growth  in  nearly  two- 
dimensional  flow  with  pressure  gradients.  The  following  development  of  Moses’s  relationships 
can  be  found  in  Hill  and  Peterson  (1992:1 19-121).  Moses’s  method  begins  with  minor  changes  to 
the  momentum  integral  equation  (Eq.2-2). 

]^[uiU-u)[dy+^]{lJ-u)dy  =  +—-—  (2-5) 

^dx  dx  ^  p  p 

where:  yi  =  any  location  within  the  boundary  layer 

'll  =  local  shear  stress  at  y  =  y/ 

To  solve  Eq  2-5,  the  quantities  5*  and  6  must  be  known,  and  8*  can  be  found  using  Eq.2- 
3a  where  the  upper  hmit  of  integration  is  the  actual  boundary  layer  thickness,  8.  Solving  for  0  is 
not  straight  forward  because  local  wall  shear  stress  can  not  be  related  to  local  pressure  gradient  in 
turbulent  flow  (Hill  and  Peterson,  1992: 119).  Instead,  Moses  integrated  Eq.2-5  from  y  =  0  to  y  = 
0.38  making  use  of  the  velocity  profiles  in  two  regions  of  the  boundary  layer:  the  nearly  laminar 
inner  layer  and  the  turbulent  outer  layer. 

A  velocity  profile  determined  by  Moses  for  each  region  was  used  to  define  a  momentum 
integral  containing  terms  related  to  the  skin  friction,  local  Reynolds  number,  freestream  velocity, 
and  kinematic  viscosity.  Each  term  in  the  equations  had  a  different  coefficient  assigned  to  it  for  a 
total  of  eight  coefficients.  These  coefficients  are  functions  that  have  to  be  calculated  separately 
during  an  iteration.  If  a  full  numerical  integration  is  not  desired,  8*  and  0  can  be  approximated  at 
discrete  points  along  a  nozzle  wall,  where  the  boundary  layer  thickness  is  known,  with  Eqs.  2-6 
and  2-7. 
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and:  Xo  =  local  wall  shear  stress 

Ap  =  pressure  drop  in  length  L 
R  =  radius  of  channel 
V  =  local  viscosity 
5  =  boundary  layer  thickness 


(2-7) 


The  boundary  layer  is  often  separated  into  three  regions  when  considering  boundary  layer 
velocity  profiles.  There  is  an  inner  layer  where  viscous  shear  dominates,  an  outer  layer  where 
turbulent  shear  dominates,  and  an  overlap  layer  where  a  smooth  transition  between  the  inner  and 
outer  layer  occurs  (White, 1991:411).  The  velocity  profile  in  a  turbulent  boundary  layer  varies 
depending  on  the  type  of  pressure  gradient  present,  but  the  log-law  using  inner  layer  variables  is 
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used  to  relate  any  type  of  pressure  gradient  to  a  single  velocity  profile.  The  log-law  with  inner 
variables  is  given  in  Eq.2-8. 


where:  v*  = 


V  H  y 


“  1  .  r. 
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and:  v  =  local  viscosity 
T  =  wall  shear  stress 
y  =  distance  from  wall 
u  =  flow  velocity  inside  boundary  layer 
k  =  0.41 


B  =  5.0 


(2-8) 


Coles  made  the  log  law  better  by  noticing  that  excess  velocity  in  the  outer  layer,  which  the 
log-law  only  adequately  predicts,  has  a  wavelike  shape  (White,1991:417).  Therefore,  a  wake 
function  was  added  to  the  log-law,  which  makes  its  velocity  profile  predictions  more  accurate. 
Defining  =  u/v*  and  y*'  =  yvVv  ,  the  log-law  with  the  addition  of  Cole’s  Law  of  the  Wake  is 
given  in  Eq.2-9. 


1  2n 

«"■  =  — ln(y'')+5  + - f\ 

K  K 


vsj 


(2-9) 


where:  11  =  — 

2 

and:  A  =  outer- variable  constant  (1.0  for  a  strong  favorable  pressure  gradient) 
B=  inner-variable  constant  (5.0) 
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n  =  Cole’s  wake  parameter 


Equation  2-9  is  complete  and  accurate  for  any  two-dimensional,  turbulent  boundary  layer 
profile,  and  if  it  is  integrated  over  the  boundary  layer,  5*  and  0  can  be  determined  from  Eq.2-10. 
Notice  that  in  this  case  the  skin-friction  coefficient  is  related  to  FI  and  X.  This  makes  the 


calculation  of  Cf  much  simpler  and  could  be  applied  in  this  form  to  other  analyses. 

5*  i  +  n 

6 


^_^_2-h3.2n-n.5n' 
6  ~  5 


(2- 10a) 

(2- 10b) 


where:  X  = 


) 


(  .  A 


=-id 

K 


^Re  ^ 


U-B+- 


in 


By  combining  the  two  methods  developed  above,  an  approximate  analysis  of  the  boundary 
layer  thickness  at  the  exit  can  be  accomplished,  but  this  is  accurate  for  a  non  blowing  wall  only. 
The  boundary  layer  would  be  thicker  on  a  blowing  wall  because  of  the  rough  porous  plate. 

k  Re 

Assuming  that  the  porous  plate  is  fully  rough  and  —  =  — Cf  and  Re^  can  be  corrected  for 

X  1000 

roughness  with  Eqs.  2-11  and  2-12  and  the  above  analysis  reaccomplished. 


Re,  =  113(1+ 03k^)e^' 
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(2-12) 


Z  =  K^ 

where:  Re^(^/;c) 

^  "  I 

and:  k  =  roughness  height 


■ 

2.87  + 1.58  log,o 

jJ 

Another  useful  analysis  is  estimating  the  boundary  layer  thickness  at  the  throat  of  a  nozzle 
because  the  thickness  is  often  too  small  to  measure  visually.  The  discharge  coefficient  is  a  well 
known  quantity  in  rocket  nozzle,  and  is  given  in  Eq.  2-13  where  R  is  defined  by  Eq  2-14. 


C^=l- 
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-2.128 
.  Y  +1 
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(2-13) 


(2-14) 


where:  y  =  ratio  of  specific  heats  (Cp/Cv) 

Re*  =  Reynolds  number  at  the  throat 
R*  =  throat  radius 

Rc  =  radius  of  curvature  of  the  throat 


It  can  be  shown  that  for  cases  where  5,  /R*«l,  the  boundary  layer  thickness  at  the  throat 
can  be  estimated  by  Eq.2-15.  Using  this  method  for  the  nozzle  used  in  this  experiment,  5,  was 
estimated  to  be  0.018  mm  (0.0007")  (See  Appendix  E).  In  this  study,  boundary  layer  thickness 
was  assumed  to  be  0.018  mm  when  the  boundary  layer  thickness  was  not  visually  measurable. 
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(2-15) 


5,  1-Cz> 

R*  2 


2.2  Transpired  Turbulent  Boundary  Layer 

One  approach  to  the  transpired  boundary  layer  was  presented  by  Kays  and  Crawford 
(1980:180-181).  They  developed  an  algebraic  Couette  flow  solution  assuming  flat  plate 
conditions.  They  found  that  the  skin  friction  at  a  given  blowing  ratio  can  be  successfully 
approximated  Eq.  2-16,  where  the  blowing  ratio  is  defined  in  Eq.  2-17. 


C^/2  ln(l-(-5^) 


"I'/m, 

C,/2 


where:  m"  =  mass  flow  rate  through  the  porous  wall 
nioo  =  freestream  mass  flow  rate 


(2-16) 

(2-17) 


Equations  2-16  and  2-17  are  solved  simultaneously  to  calculate  the  change  in  C//  2  with 
blowing.  Once  the  new  C//  2  is  known,  roughness  can  be  taken  into  account,  and  other  methods 
to  approximate  the  boundary  layer  thickness  with  transpiration  can  be  used. 

Another  approach  to  determining  the  effect  of  transpiration  on  the  boundary  layer  is  to 
modify  the  log-law  with  another  logarithmic  term,  which  was  accomplished  by  Stevenson 
(White,  199 1:435).  His  modified  log-law,  defining  Vw^  =  Vw/v*,  is  given  in  Eq.2-18.  With  this 
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equation,  law  of  the  wall  plots  can  be  created  that  include  blowing  and  allow  an  estimate  of  Q 
reduction  with  blowing  to  be  obtained.  After  correcting  Q,  the  boundary  layer  thickness  with 


blowing  can  be  approximated  using  the  previous  analyses  by  Moses  and  Coles. 

^[(l+v,V)^-l]--iln(/)  +  B 

V^,  1C 


(2-18) 


It  should  be  noted  that  porous  plates  cause  roughness  and  disturbance  effects  on  the 
boundary  layer,  but  it  was  assumed  that  roughness  caused  the  same  effect  in  each  case,  and  any 
additional  growth  beyond  that  was  a  result  of  blowing. 
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III.  Experimental  Setup 


3.1  Shock  Tunnel  Configuration 

According  to  shock  wave  theory,  when  a  shock  wave  passes  through  an  undisturbed  fluid 
it  increases  the  pressure  and  temperature  of  the  fluid  behind  it  (Shapiro,1954:116)  as  well  as 
accelerating  the  fluid  in  the  direction  of  the  shock  wave  (Shapiro,1954:998).  If  the  shock  wave 
reflects  with  equal  strength  from  a  straight  wall  the  same  effect  occurs  but  the  acceleration  is  in 
the  opposite  direction  causing  the  fluid’s  net  velocity  to  be  zero  behind  the  shock  with  even  higher 
pressure  and  temperature.  This  is  the  reason  shock  tunnel  facilities  are  suited  for  this  kind  of 
experiment  (Shapiro,  1954:1021,1008).  The  “pocket”  of  zero  velocity,  high  pressure,  high 
temperature  air  resulting  from  the  reflected  shock  provided  the  simulated  chamber  conditions  for 
the  present  study. 

The  AFIT  low  pressure  shock  tunnel  was  the  faeility  used  for  this  experiment,  and  Figure 
3.1  shows  a  schematic  of  the  facihty  including  all  external  connections.  The  shock  tunnel  has  a 
1.22  m  (4')  driver  section,  a  4.88m  (16')  driven  section,  a  1.22  m  (4')  test  section,  and  is  20.3  cm 
(8")  tall  by  10.2  cm  (4")  wide.  This  experiment  used  a  0.18mm  (0.007")  mylar  diaphragm  to 
separate  the  high  pressure  section  from  the  low  pressure  section. 

The  tunnel  was  operated  with  an  air/helium  driver  mixture,  which  increased  the  strength  of 
the  shock  and  the  chamber  conditions  to  a  level  higher  than  an  air  only  driver  would  have 
produced  (Shapiro,  1954:9 18).  Using  an  air/He  driver  mixture  changes  the  ratio  of  specific  heats 
(y)  in  the  driver,  and  this  value  must  be  known  to  accurately  calculate  the  chamber  conditions 
after  shock  reflection.  Appendix  A  shows  an  example  calculation  of  y  for  the  air/He  mixture. 
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Figure  3.1  -  Shock  Tunnel  Configuration 


In  this  experiment,  simulated  chamber  conditions  were  produced  when  the  shock  wave 
reflected  from  the  converging  section  of  a  Mach  2.0  characteristic  nozzle  located  in  the  test 
section  so  that  the  driven  section  was  5.49  m  (18')  long.  Knowing  the  driven  section  length, 
pressure  and  temperature  as  well  as  the  driver  pressure,  temperature,  and  y,  the  chamber 
conditions  could  be  predicted  according  to  shock  tunnel  theory  (Shapiro,  1954: 1008-1009). 
Program  SHOCKTUN.FOR  (Bowersox,1990)  calculated  the  conditions  in  each  region  of  the 
shock  tunnel  for  given  test  conditions,  and  although  test  conditions  varied  shghtly  from  run  to 
run.  Table  3.1  shows  the  design  driver  conditions,  analytical  chamber  conditions,  and  analytical 
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run  time  for  the  experiment.  Appendix  B  contains  instructions  for  using  SHOCKTUN.FOR.  It 
was  expected  that  the  actual  chamber  conditions  would  be  less  than  the  values  shown  in  Table  3.1 
because  the  shock  wave  was  not  perfectly  reflected  from  the  converging  walls  of  nozzle  and  some 
of  the  shock  wave  was  swallowed  by  the  nozzle. 


Table  3.1  -  Summary  of  Expected  Test  Conditions 


Condition 

Driver  Pressure 

0.58  MPa  (84.4  psi) 

Air 

0.10  MPa  (14.3  psi) 

Helium 

0.48  MPa  (70  psi) 

Driven  Pressure 

0.10  MPa  (14.3  psi) 

P5 

0.71  MPa  (103.1  psi) 

T5 

546.0  K 

Run  Time 

8.28  ms 

The  only  other  configuration  note  on  the  shock  tunnel  is  that  the  downstream  end  of  the 
tunnel  was  left  open  to  prevent  undesirable  reflections.  Therefore,  shock  waves  were  allowed  to 
pass  out  of  the  shock  tunnel  and  into  the  room.  This  was  acceptable  because  the  exit  pressure  of 
the  nozzle  was  close  enough  to  atmospheric  pressure  that  flow  separation  due  to  overexpansion 
and  shock  waves  in  the  nozzle  due  to  a  deficient  pressure  ratio  were  avoided. 


3.2  Test  article 

3.2.1  General  Description 

The  test  article  used  for  this  experiment  was  a  Mach  2.0  characteristic  nozzle  designed 
using  three  computer  programs  written  for  nozzle  design  by  Dr.  Rodney  Bowersox  (1990). 
Figure  3.2  shows  the  features  and  dimensions  of  the  model. 
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(b)  Non  Blowing  Side  View 


(c)  Blowing  Side  View 
Figure  3.2  -  Nozzle  Design  w/  Dimensions 
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Some  nozzle  features  were  chosen  arbitrarily  because  the  model  only  needed  to  be  a 
supersonic  nozzle  appropriate  for  a  shock  tunnel  study  of  boundary  layer  thickness  and  heat 
transfer  reduction  in  a  transpiration  cooled  nozzle.  Mach  2.0  was  chosen  because  it  would  clearly 
be  supersonic  even  with  losses  included,  the  pressure  ratio  for  choked  flow  could  be  achieved 
without  drawing  a  vacuum  in  the  shock  tunnel,  and  a  diffuser  was  not  necessary  to  avoid  flow 
separation  and  shock  waves  in  the  nozzle. 

The  throat  size  was  chosen  based  on  the  requirement  that  it  be  small  compared  to  the  size 
of  the  shock  tunnel  so  that  the  converging  section  of  the  nozzle  would  approximate  a  straight  wall 
for  the  shock  wave  to  reflect  from.  Zuppan  suggested  in  his  work  that  a  throat/shock  tunnel 
height  ratio  of  10.65  to  13.3  was  appropriate  for  good  shock  wave  reflection,  and  in  his  work  he 
used  a  throat  height  of  1.91  cm  (0.75")  (Zuppan,  1965:Preface).  Because  the  same  shock  tunnel 
was  being  used  for  this  experiment,  Zuppan’ s  choice  seemed  appropriate  and  was  adopted  here. 
In  addition,  the  model  was  designed  to  be  as  wide  as  the  shock  tunnel  (4"),  which  made  the  flow 
two-dimensional. 

One  side  of  the  nozzle  was  designed  for  transpiration  cooling  while  the  other  remained 
impermeable  to  act  as  a  control  in  the  experiment.  Instruments  were  inserted  along  the  non 
blowing  side  of  the  nozzle  through  a  circular  cavity  drilled  into  it  from  the  bottom  (see  Fig  3.2b). 
Also,  chamber  pressure  was  measured  from  the  non  blowing  side  and  a  deep  channel  had  to  be 
drilled  to  insert  this  transducer. 

The  cavity  on  the  blowing  side  was  square  and  went  all  the  way  through  the  waU.  During 
blowing  this  cavity  acted  as  a  constant  pressure  plenum  to  ensure  that  the  porous  material 
experienced  the  same  pressure  throughout  the  test  as  well  as  providing  space  to  insert  instmments 
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along  the  wall.  From  a  point  1.27  cm  (0.5")  before  the  throat  to  the  nozzle  exit,  the  wall  was 


recessed  by  0.635  cm  (0.25")  creating  a  contoured  shelf  on  the  blowing  side,  and  the  four  layers 
of  the  selected  porous  material  (see  Appendix  F)  were  installed.  Each  layer  of  porous  material 
was  also  shaped  to  the  nozzle  contour  individually  to  minimize  the  unevenness  at  the  porous 
material/aluminum  seam,  which  could  trip  the  flow  and  disrupt  the  flow  field  downstream  (see 
Figure  3.2c).  Even  with  this  precaution  there  was  still  enough  unevenness  at  that  point  to  trip  the 
flow  (see  Figure  3.3).  Because  the  flow  tripped  in  the  converging  section  where  the  flow  was 
subsonic,  it  quickly  reattached  to  the  wall. 

The  selected  porous  material  was  1.57  mm  (0.062")  thick  316L  sintered  stainless  steel 
with  2  micron  porosity,  which  was  selected  based  on  strength,  mass  flow  rate,  and  pressure  drop 
considerations.  The  porous  surface  was  clearly  more  rough  than  the  non  blowing  side,  and  the 
boundary  layer  was  expected  to  be  thicker  on  the  blowing  side  in  all  cases. 


Figure  3.3  -  Tripped  Flow  at  Porous  Material/Aluminum  Seam 
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3.2.2  Nozzle  Design  Tools 


As  mentioned,  three  computer  programs  were  used  to  design  the  Mach  2.0  nozzle,  and  the 
mechanics  of  each  program  used  will  be  treated  briefly.  Other  specifics  such  as  input  and  output 
information  on  each  program  can  be  found  in  Appendix  B. 

Program  NOZ2.FOR  (Bowersox,1990)  constructed  the  contoured  and  minimum  length 
diverging  section  of  the  Mach  2.0  nozzle  shape  using  the  method  of  characteristics  with  two 
characteristic  reflections.  Program  BLCORR.FOR  (Bowersox,1990)  calculated  the  analytical 
boundary  layer  thickness  at  the  nozzle  exit  and  adjusted  the  entire  diverging  section  contour  to 
compensate  while  maintaining  the  desired  Mach  number.  These  adjustments  change  the  throat 
and  exit  geometries  slightly  resulting  in  a  nozzle  that  was  not  ideal  according  to  the  isentropic 
flow  tables  in  Shapiro  but  was  very  close.  Figure  3.2  shows  the  nozzle  after  adjustments  made  by 
this  program,  and  the  actual  design  had  A/A*  =  1.738  and  Me  =  2.03. 

Program  WALL.FOR  (Bowersox,1990)  constructed  the  converging  section  of  the  nozzle 
and  attached  it  to  the  diverging  section.  This  program  was  used  to  design  the  converging  section 
of  the  nozzle  in  such  a  way  as  to  enhance  the  shock  reflection  (i.e.,  steep  walls),  yet  maintain  a 
reasonable  convergence  rate.  There  is  a  parameter  in  the  program  that  defines  the  length  of  the 
sloped  part  of  the  converging  section,  and  it  was  chosen  to  be  small  to  produce  steep  converging 
walls. 

3.3  Tnstmmentation  and  Data  Collection 

Ambient  pressure  was  measured  with  a  lab  barometer,  and  thermocouples  were  used  to 
measure  ambient  and  plenum  temperatures.  Other  pressure  measurements  were  collected  using 
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10  Endevco  Model  8530A-100  transducers  (range:  0-100  psia)  and  one  Endevco  Model  8514A- 
50  transducer  (range:  0-50  psia).  Their  signals  were  amplified  using  Endevco  Model  4423  Signal 
Conditioners,  and  transducer  calibration  was  accomplished  using  a  standard  calibration  setup. 
Calibration  data  is  included  in  Appendix  C.  Two  Model  8530  transducers  were  plugged  into  the 
top  of  the  shock  tunnel  itself  for  shadowgraph  timing  purposes  (henceforth  referred  to  tube 
transducer  1,  which  is  further  upstream  of  the  nozzle,  and  tube  transducer  2).  The  time  between 
the  pressure  jumps  caused  by  the  shock  wave  passing  by  the  transducers  was  the  reference  time 
for  triggering  the  light  source  (see  Section  3.4).  One  Model  8530  transducer  was  placed  in  the 
converging  section  of  the  model  to  measure  chamber  pressure,  one  was  placed  in  the  plenum  to 
measure  the  pressure  of  injected  air,  and  the  other  six  Model  8530s  transducers  were  placed  along 
the  nozzle  in  pairs  directly  opposite  each  other  to  measure  static  pressure  along  the  wall.  The 
Model  8514  transducer  was  mounted  inside  of  a  pitot  probe  and  used  to  measure  exit  pressure. 
All  Signal  Conditioners  were  set  to  a  gain  of  20X  except  those  amplifying  the  two  tube  transducer 
signals,  which  were  set  to  a  gain  of  50X.  Figures  3.2b  and  3.2c  sununarize  the  transducer 
placement  along  the  nozzle. 

The  Nicolet  500  Data  Collection  system  is  a  commercial  data  collection  system  that 
operates  in  the  Windows  environment.  The  equipment  and  documentation  was  located  in  the  lab 
area.  Pressure  data  taken  with  the  Nicolet  was  sampled  at  500  kHz 

3.4  Shadowgraph  Flow  Visualization  System 

A  shadowgraph  is  a  photograph,  which  can  be  taken  for  flow  visualization.  A 
shadowgraph  system  requires  three  basic  items:  an  aperture,  a  light  source,  and  a  mirror.  In 
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general,  the  mirror  is  placed  so  that  light  reflected  from  it  will  Uluminate  the  area  to  be 
photographed.  The  light  source,  the  mirror,  and  the  desired  image  should  be  at  the  same  height  to 
ensure  that  the  light  passes  perpendicularly  through  the  flow.  The  light  source  is  placed  at  the 
mirror’s  focal  length  from  the  mirror  so  that  when  the  light  is  reflected  from  the  mirror  it  will  be  a 
uniform  column.  A  shallow  angle  between  the  light  source  and  reflected  column  helps  to  ensure 
that  the  column  of  light  is  uniform.  The  column  passes  through  the  flow  being  photographed  and 
exposes  the  film,  and  because  of  density  gradients  in  the  flow  field,  features  such  as  shock  waves 
and  boundary  layers  are  visible.  Figure  3.4  shows  the  basic  setup  used  in  this  experiment. 


For  tests  with  long  run  times,  this  system  can  be  set  up  and  the  film  exposed  by  hand  at 
any  point  during  the  test.  Run  times  in  a  shock  tunnel,  however,  are  on  the  order  of  milliseconds, 
and  it  is  impossible  to  trigger  the  light  source  by  hand  during  this  time  especially  when 
shadowgraphs  at  specific  points  during  the  run  time  are  desired.  Therefore,  an  electrical  system 
was  employed  to  trigger  the  light  source  at  a  preset  time  measured  in  microseconds  after  the 
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shock  wave  passed  tube  transducer  2.  The  shadowgraph  system  complete  with  a  timing  circuit 
was  setup  to  capture  events  in  the  nozzle  before,  during,  and  after  its  run  time,  and  the  system 
was  demonstrated  to  be  very  reliable  and  accurate  during  this  experiment. 


Figure  3.5  -  Pitot  Probe  Bent  Out  of  Flow 


For  example,  at  one  point  the  pitot  probe  was  not  collecting  usable  data  when  it  should 
have  been.  A  shadowgraph  was  used  to  capture  the  pitot  probe  bent  out  of  the  flow  (see  Fig.  3.5) 
and  hence,  reading  only  ambient  pressure.  This  discovery  resulted  in  stiffening  modifications  to 
the  pitot  probe  that  made  all  pitot  measurements  more  reliable  as  well  as  saving  a  lot  of  time  and 
frustration  in  determining  the  source  of  the  poor  data.  Table  3.2  summarizes  the  electrical 
components  of  the  system  and  their  functions,  and  Figure  3.6  shows  the  system  in  block  diagram 
form. 
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Table  3.2  -  Timing  Circuit  Components  w/  function 


Figure  3.6  -  Timing  Circuit  Block  Diagram 


The  spark  lamp  trigger  time  delay  began  when  the  shock  wave  passed  by  tube  transducer 
2,  and  the  total  time  delay  is  based  on  the  measured  amount  of  time  (measured  by  the  Model  436) 
between  tunnel  transducers  1  and  2.  From  Figure  3.6,  then,  it  can  be  seen  that  the  total  delay 
period  until  the  light  source  is  triggered  can  be  set  to  any  time  by  satisfying  Eq.  3-1. 
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TD^  —  +  ^^453 


(3-1) 


where:  TD,  =  total  time  delay  (^s) 

TD436  =  measured  time  delay  on  Model  436  (p,s) 

TD453  =  time  delay  set  on  Model  453  ()is) 

In  practice,  TD436  was  measured,  and  7D,  was  the  target  total  time  delay  based  on  TD436 
and  the  event  to  be  captured.  Therefore,  the  value  really  solved  for  in  Eq.3-1  was  TD453.  Using 
this  system  required  at  least  one  run  at  the  desired  test  condition  without  a  photograph  but  with 
the  system  running  otherwise  to  measure  TD436-  It  was  known  that  the  distance  between  the  two 
tube  transducers  was  0.71  m  (28"),  and  therefore,  once  TD436  was  known  and  assuming  the  shock 
wave  propagation  speed  was  constant,  Tdt  could  be  determined  by  Eq.3-2. 


7^^436  _  TP, 

28  (28-i-x) 

where:  x  =  distance  past  the  second  tunnel  transducer  (in) 


(3-2) 


The  proportional  relationship  of  Eq.3-2  was  used  by  adjusting  TD453  appropriately  and 
was  probably  most  useful  as  a  tool  for  capturing  an  initial  event.  In  this  case,  x  was  chosen  to  be 
some  point  along  the  model.  The  assumption  that  the  shock  wave  was  moving  at  a  constant 
speed  held  up  for  test  runs  with  an  air  only  driver  as  shock  waves  were  captured  at  various 
locations  along  the  nozzle  with  great  accuracy.  The  linear  relationship  did  not  hold  up  for  the 
air/He  driver,  which  indicated  that  the  shock  wave  had  actually  slowed  down  slightly  after  passing 
by  the  tube  transducers.  A  correction  of  250  |i.s  was  added  to  TD453  to  capture  the  initial  event. 
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Once  one  event  was  captured,  TD453  was  simply  incremented  forward  and  backward  until  the 
desired  range  of  events  was  captured.  This  range  might  include  such  highlights  as  nozzle  start, 
established  flow,  and  unstart.  Test  conditions  changed  from  run  to  run,  which  affected  TD436,  and 
changed  the  delay  slightly  for  each  run.  Experience  in  this  experiment  indicated  that  if  the 
conditions  did  not  vary  greatly  (i.e.,  driver  pressure  within  1  in  Hg)  then  TD436  varied  by  10  |xs  at 
most.  Because  the  run  times  were  on  the  order  of  milliseconds,  this  difference  never  caused  the 
desired  event  to  be  missed. 

The  general  procedure  for  using  the  timing  system,  then,  was  to  run  a  test  to  measure  the 
approximate  TD436,  set  TD4S3  to  an  appropriate  value  for  capturing  the  desired  event,  run  the  test, 
and  reset  the  system.  The  most  important  component  to  reset  was  the  homemade  gray  box 
because  it  was  a  high  voltage  component  that  stored  its  -250V  pulse  in  a  capacitor  that  would 
quickly  bum  up  if  not  reset.  The  box  was  reset  by  pushing  the  small  orange  button,  and  the 
orange  light  should  come  on.  Shadowgraph  procedures  were  incorporated  into  the  shock  tunnel 
operation  procedures  because  the  two  procedures  did  occur  simultaneously  (see  Appendix  D  for 
procedures). 

Documentation  is  available  on  all  components  with  the  exception  of  the  homemade  gray 
box  and  divider  circuit.  The  homemade  gray  box  has  no  documentation  except  what  is  provided 
here,  and  Figure  3.7  shows  a  schematic  of  the  divider  circuit. 
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Figure  3.7  -  Schematic  of  Voltage  Divider  Circuit 


3.5  Procedural  Notes 

There  are  a  few  miscellaneous  items  that  need  to  be  discussed  on  the  setup  and  operation 
of  the  AFIT  low  pressure  shock  tunnel.  First,  the  driver  pressure  is  measured  in  in  Hg  when 
mnning  a  test  and  unit  conversions  should  be  double  checked  for  accuracy  to  ensure  the  correct 
test  conditions  are  being  used.  Also,  the  diaphragm  stretches  as  the  driver  is  pressurized,  which  is 
expected,  but  this  stretching  can  cause  the  shifting  of  the  diaphragm  if  it  is  not  secured  tightly  into 
place.  This  condition  tended  to  cause  premature  bursting  of  the  diaphragm.  Periodic  clicking 
sounds  as  the  driver  is  pressurized  is  a  warning  signal  for  this  condition. 

Next,  it  was  found  that  after  a  test  using  an  air/He  driver  the  entire  tunnel  was  filled  with 
the  air/He  mixture,  and  if  a  test  was  run  immediately  after  the  previous  one  this  new  mixture  in 
the  driven  section  affected  the  shock  speed.  Therefore,  a  five  minute  helium  purge  period  was 
instituted  between  each  test  in  which  house  air  is  forced  through  the  tunnel  to  evacuate  excess 
helium.  Shock  speed  results  were  very  consistent  after  the  procedure  was  implemented. 
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Also,  when  taking  shadowgraphs,  it  is  important  to  remember  that  the  film  is  exposed 
throughout  the  test  so  the  test  area  must  be  very  dark  to  avoid  unwanted  or  premature  exposure. 

Lastly,  if  the  shock  tunnel  is  run  with  the  end  open  as  in  this  experiment,  ear  protection 
should  be  worn  because  it  is  very  loud.  Also,  it  is  best  if  the  lab  door  is  closed,  and  all  personnel 
in  the  room  should  be  warned  when  a  test  is  about  to  be  run  before  every  test. 
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IV.  Results  and  Discussion 


4.1  Baseline  Results 

4.1.1  Actual  Test  Conditions 


Recall  that  it  was  expected  that  the  actual  chamber  conditions  would  be  less  than 
SHOCKTUN.FOR  predictions  because  the  shock  wave  was  reflecting  from  a  converging  nozzle 
section  instead  of  a  straight  wall.  In  fact,  chamber  pressure  and  temperature  were  32%  and  12% 
less  than  predictions  respectively.  Temperature  inside  the  shock  tunnel  was  not  actually 
measured,  however,  but  was  estimated  with  SHOCKTUN.FOR  according  to  what  the  chamber 
pressure  would  have  been  if  the  measured  Po  were  the  actual  chamber  condition  after  a  perfect 
reflection.  Chamber  temperature  (To)  was  estimated  at  other  chamber  conditions  using  the 
proportional  relationship. 


To  =  (0.8699^ 


p  ^ 
M) 

103.1 


%06716)- (545.978  K) , 


which  was  based  on  the  difference  between  predicted  and  actual  conditions.  This  was  acceptable 
because  chamber  pressure  did  not  vary  outside  the  linear  range.  Table  4.1  sununarizes  the  actual 
test  conditions  for  the  experiment. 


Table  4.1  -  Summary  of  Actual  Test  Cond 

itions 

Condition 

%  Error 

Driver  Pressure 

580  kPa  (84.3  psi) 

— 

Air 

98.6  kPa(14.3  psi) 

— 

Helium 

482.5  kPa  (70  psi) 

— 

Driven  Pressure 

98.6  kPa(14.3  psi) 

___ 

P5 

482.3  kPa  (69.983  psi) 

±1.04% 

T5 

480.1  K 

±1.04% 

Shock  Speed 

527.54  m/s 

±0.11% 

Run  Time 

1.938  ms 

±3.92% 

Notice  from  the  error  column  of  Table  4.1  that  conditions  were  very  repeatable,  especially 
the  shock  strength,  which  was  proportional  to  the  shock  speed.  The  shock  speed  was  calculated 
based  on  the  time  it  took  the  shock  wave  to  pass  the  tunnel  transducers.  Test  conditions  were 
consistent  enough  that  the  data  was  reduced  without  normalization  in  all  cases  except  the  specific 
impulse  analysis. 

The  run  time  in  Table  4.1  was  76%  lower  than  predicted  (See  Table  3.1).  This  was  due, 
in  part,  to  losses,  but  was  mostly  due  to  the  different  definitions  of  run  time  used  in  the  two 
tables.  In  Table  3.1,  run  time  was  defined  as  the  time  between  shock  wave  reflection  and  the 
contact  discontinuity  reaching  the  end  wall.  In  Table  4.1,  the  run  time  was  defined  as  the  time 
after  shock  wave  reflection  that  the  nozzle  maintained  steady  state  flow  conditions.  The  end  of 
the  run  time,  then,  was  deflned  as  the  time  when  the  nozzle  pressure  ratio  became  deficient 
because  of  chamber  pressure  depletion  and  unstart  occurred. 

Shadowgraphs  were  vital  to  bracketing  the  actual  run  time  of  the  nozzle  because  they 
showed  very  clearly  the  start  and  unstart  of  the  nozzle.  Because  the  timing  of  each  photograph 
was  known  (i.e.,  the  timing  circuit.  Sec  3.4),  the  run  time  could  be  estimated.  Run  times  were 
measured  accurately  from  pitot  probe  data  that  also  showed  clearly  the  beginning  and  end  of 
steady  state  flow  (See  Figure  4.  If). 

4. 1 .2  Characterization  of  the  Flow  Field 

Shadowgraphs  were  taken  at  several  times  during  the  nozzle  starting  process  to  visualize 
and  understand  the  events  leading  to  established  flow  in  the  nozzle.  The  photographs  did  point 
out  that  the  flow  established  at  each  point  along  the  nozzle  about  1  ms  after  the  shock  wave 
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Figure  4.1  -  Characterization  of  Flow  Dynamics 


passed  by,  and  this  was  a  data  landmark.  Figure  4.1  is  a  pictorial  summary  of  a  non  blowing  test 
run  from  initial  shock  wave  reflection  to  unstart,  which  took  place  in  about  3.5  ms. 

In  the  photos,  the  flow  is  propagating  from  left  to  right,  and  the  upper  wall  is  the  blowing 
side  of  the  nozzle.  This  photograph  orientation  is  maintained  throughout  this  report.  The 
pressure  traces  shown  are  centerline  pitot  probe  and  chamber  pressure,  and  they  are 
representative  of  runs  throughout  the  experiment  (See  Figure  4.  If).  Constant  chamber  pressure 
over  all  test  runs  was  attributed  to  test  condition  consistency.  All  pitot  traces  showed  the  same 
general  form  except  that  run  times  were  shorter  when  the  probe  was  close  to  the  wall  because  of 
the  unstart  geometry  (See  Figure  4.1e).  These  shadowgraphs  are  also  typical  because  blowing 
had  no  visible  effect  on  the  start  up  dynamics.  Figure  4. 1  also  helps  to  characterize  the  difference 
between  the  non  blowing  side  and  the  blowing  side  under  baseline  conditions. 

4.1.3  Non  Blowing  vs  Blowing  Side 

Figure  4.1c  shows  the  entire  nozzle  at  steady  state  conditions  including  boundary  layer 
establishment  and  growth  along  the  nozzle.  In  the  photo,  the  thin,  bright  region  near  each  wall  is 
the  boundary  layer  for  that  side  of  the  nozzle,  and  there  is  clearly  a  visible  difference  in  thickness 
between  the  two  sides.  Using  a  caliper  micrometer,  boundary  layer  thickness  measurements  were 
made  from  the  actual  photographs  and  from  projected  images  of  the  photographs.  Several  trial 
measurements  were  taken,  and  the  statistical  average  of  those  measurements  for  each  case  was 
defined  as  the  boundary  layer  thickness.  Measurements  were  made  the  same  way  at  positions 
along  the  nozzle  to  show  the  rale  jaf  boundary  layer  growth.  The  exit  boundary  layer  thickness 
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was  of  the  most  interest  because  this  thickness  decreased  the  effective  exit  area  of  the  nozzle, 
which  affected  its  performance. 

For  the  baseline  case,  the  exit  boundary  layer  thickness  was  213%  greater  on  the  blowing 
side  than  on  the  non  blowing  side.  Table  4.2  compares  the  exit  boundary  layer  thicknesses  (He 
represents  exit  height),  and  Figure  4.2  compares  the  boundary  layer  thickness  on  both  sides  as  a 
function  of  distance  from  the  throat. 


Table  4*2  -  Comparing  Non  Blowing  and  Blowing  Sides  of  Nozzle 
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Clearly,  the  mere  presence  of  a  rough  porous  wall  had  an  effect  on  the  boundary  layer,  but 
this  was  not  a  true  baseline  case.  Evidence  that  the  roughness  of  the  porous  material  was  the 
major  cause  of  boundary  layer  growth  in  the  baseline  case  is  provided  by  calculating  k"^.  This 
value  is  a  Reynolds  number  based  on  the  average  roughness  of  the  surface.  Taking  the  average 
roughness  to  be  on  the  order  of  the  porosity  of  the  material  (2  micron  porosity),  k"^  was  calculated 
to  he  76  at  the  throat  and  58  at  the  exit.  When  k^  is  greater  than  60,  the  surface  is  considered  to 
be  fully  rough,  which  means  that,  in  this  case,  the  surface  was  fully  rough  over  almost  the  entire 
porous  surface.  Based  on  this  evidence,  roughness  of  the  porous  material  was  the  major 
contributor  to  boundary  layer  growth  in  the  baseline  case.  It  must  be  stated  here,  however,  that  a 
true  baseline  case  was  not  observed  in  the  experiment,  and  this  also  had  an  effect  on  the  baseline 
results. 

When  there  was  not  active  injection  through  the  porous  material,  the  plenum  pressure  was 
equal  to  atmospheric  pressure.  The  pressure  gradient  still  existed  in  the  nozzle  during  a  run  so 
pressure  drops  still  occurred  through  the  porous  material.  Therefore,  even  without  active 
injection,  there  was  a  very  small  amount  of  suction  near  the  throat  of  the  nozzle  where  the  wall 
pressure  was  greater  than  atmospheric  and  a  very  small  amount  of  blowing  near  the  end  of  nozzle 
(B=.000128±7.2%)  as  the  wall  pressure  dipped  below  atmospheric.  Both  of  these  effects  could 
also  have  contributed  to  the  boundary  layer  growth  on  the  blowing  side. 

The  presence  of  suction  and  blowing  in  the  baseline  case  with  no  active  injection  imphes 
that  blowing  was  always  greatest  near  the  end  of  the  nozzle  where  the  pressure  drop  was  greatest 
and  cooling  was  needed  least.  Lack  of  flow  control  through  porous  material  is  a  major  criticism 
of  transpiration  cooling  (May  and  Burkhardt,  1991:1),  and  it  certainly  was  illustrated  here  because 
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maximum  blowing  near  the  exit  affected  the  boundary  layer  thickness  results  of  this  experiment. 
While  these  measurements  and  observations  point  out  the  differences  between  the  two  sides  of  the 
nozzle,  they  also  point  out  that  even  under  baseline  conditions,  the  boundary  layer  over  a  porous 
plate  is  a  complex  phenomenon. 

In  a  rocket  nozzle  with  porous  nozzle  walls,  the  boundary  layer  would  be  thicker  than  if 
the  walls  were  impermeable  as  shown  here.  Therefore,  additional  boundary  layer  growth  in  the 
blowing  cases  of  the  experiment  was  measured  relative  to  the  baseline  blowing  side  boundary 
layer  thickness. 

4. 1 .5  Exit  Mach  Number  Profile 

Exit  Mach  number  is  a  measure  of  performance  because  it  is  directly  related  to  exit 
velocity.  Also,  a  uniform  exit  Mach  number  profile  means  uniformly  axial  exit  flow,  which  is 
desired  for  wind  tunnels.  Exit  Mach  data  was  collected  with  a  pitot  probe  placed  at  several 
positions  across  the  exit  plane. 

The  exit  Mach  number  profile  was  calculated  by  averaging  the  total  pressure  ratio  {pta|p^^) 
in  front  of  the  pitot  probe  over  the  run  time,  and  calculating  Mach  number  using  the  gas  dynamic 
equations  (1953:623).  Table  4.3  summarizes  the  results. 


Table  4.3  -  Summary  of  Baseline  Exit  Mach  Number  Data 


Position 

Mach  Number 

Error 

Average  across  proflle 

1.886 

±1.07% 

Centerline 

1.92 

— 

Avg  in  non  blowing  region 

1.88 

±2.18% 

Avg  in  blowing  region 

1.88 

±1.56% 
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The  averages  referred  to  in  Table  4.3  and  throughout  this  report  refer  to  spatial  averages 
across  the  exit  plane,  not  the  average  over  several  runs.  Therefore,  these  parameters  were  chosen 
because  they  were,  statistically,  better  indicators  of  the  effects  of  blowing  on  exit  Mach  number 
than  a  single  data  points  near  the  walls.  The  non  blowing  and  blowing  regions  were  defined  such 
that  all  data  points  from  the  centerline  to  the  respective  wall  excluding  centerline  data  points  were 
included  in  the  calculation  of  the  average  Mach  number.  The  choice  of  these  parameter  resulted 
in  a  symmetric  distribution  across  the  exit  plane  in  the  original  baseline  case  (Table  4.3). 
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Figure  4.3  -  Baseline  Exit  Mach  Number  Profile 


As  mentioned  earlier,  the  baseline  data  did  not  represent  a  true  baseline  case,  which  would 
have  required  a  nozzle  with  two  impermeable  walls.  The  suction  and  blowing  present  in  the 
baseline  case  did  slightly  affect  the  baseline  exit  Mach  number  profile.  A  true  baseline  case  was 
created  by  assuming  that  the  Mach  number  profile  near  the  non  blowing  wall  would  be  the  same 
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near  any  impermeable  wall.  The  data  near  the  non  blowing  wall  was  mirrored  above  and  below 
the  nozzle  centerline  to  create  a  true  baseline  case  for  exit  Mach  number  profile.  Therefore,  six 
exit  Mach  number  cases  are  represented  here:  one  baseline  and  five  with  blowing.  The  blowing 
case  that  had  no  active  injection  is  referred  to  as  the  Blowing  0  case  for  the  remainder  of  this 
report.  The  new  baseline  exit  Mach  number  profile  using  the  mirrored  data  is  shown  in  Figure 
4.3. 


Figure  4.4  -  Pitot  Probe  as  Part  of  the  Wall 


Many  more  data  points  were  taken  for  the  baseline  case  than  for  the  blowing  cases 
because  a  very  complete  picture  of  the  baseline  exit  Mach  number  profile  was  desired.  Also, 
decreases  in  exit  Mach  number  with  blowing  were  expected  to  occur  most  dramatically  near  the 
walls  so  less  freestream  data  points  were  deemed  necessary.  Unfortunately,  the  pitot  probe  was 
necessarily  large  compared  to  the  exit  height  of  the  nozzle,  which  made  fine  measurements  inside 
the  boundary  layer  itself  impossible.  In  fact,  the  probe  had  to  be  placed  at  least  3.37mm  (0. 1328") 


from  the  wall  or  the  wall/probe  interaction  contaminated  the  pitot  pressure  measurements  (See 
Figure  4.4).  Therefore,  data  for  exit  Mach  number  immediately  next  to  either  wall  was  not 
acquired. 

4.2  Blowing  Cases 
4.2.1  Definition 

Four  cases  with  active  injection  were  studied,  and  they  were  defined  by  the  target  plenum 
pressure  used  in  each  case.  The  minimum  plenum  pressure  (259  kPa)  was  the  lowest  possible 
pressure  that  produced  positive  blowing  along  the  entire  porous  area.  The  maximum  plenum 
pressure  was  set  based  on  hardware  limitations. 

A  regulator  was  necessarily  connected  to  the  air  bottle  used  for  injection,  and  even  with  it 
and  the  air  bottle  completely  open,  the  regulator  was  mass  flow  limited.  Therefore,  the  plenum 
pressure  could  only  be  maintained  at  382  kPa,  which  became  the  maximum  plenum  pressure.  The 
other  plenum  pressures  were  chosen  so  that  the  intermediate  cases  would  be  evenly  spaced 
between  the  minimum  and  maximum  plenum  pressures.  The  test  plenum  conditions  were  very 
repeatable  at  each  blowing  level,  and  the  variations  in  the  plenum  conditions  were  small  enough  to 
have  very  Uttle  effect  on  the  blowing  ratios.  Table  4.4  summarizes  the  blowing  conditions.  The 
Blowing  0  case  was  included  here  because  blowing  was  present  at  atmospheric  plenum 
conditions,  which  made  it  the  limiting  case  of  blowing  as  mentioned  earlier. 

Blowing  ratios  were  calculated  by  Eq.1-1  with  program  BLOWDOT.FOR  in  which  mass 
flow  rates  were  calculated,  converted  to  the  blowing  parameter,  and  summed  at  increments  along 
the  porous  area.  Calculations  were  based  on  the  actual  x-y  coordinates  of  the  nozzle,  measured 
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static  pressure  at  each  coordinate,  and  plenum  conditions.  Each  increment  was  centered  around 
an  X  coordinate  so  that  the  summation  is  more  accurate.  Only  positive  blowing  was  considered  in 
the  blowing  ratio,  but  the  suction  mass  flow  rate  was  also  calculated  for  later  analysis.  Table  4.5 
summarizes  the  blowing  ratios  obtained  in  this  experiment. 


Table  4.4  -  Summary  of  Blowing  Conditions 


Case 

Plenum  Pressure 

Error 

Plenum  Temp 

Error 

Blowing  0 

98.5  kPa  (14.29  psi) 

±0.09% 

295.18  K 

±0.04% 

Blowing  1 

259  kPa  (37.63  psi) 

±1.19% 

291.53  K 

±1.10% 

Blowing  2 

314  kPa  (45.57  psi) 

±0.44% 

288.62  K 

±0.25% 

Blowing  3 

348  kPa  (50.55  psi) 

±3.85% 

287.61  K 

±0.35% 

Blowing  4 

382  kPa  (55.36  psi) 

±0.96% 

288.48  K 

±f).75% 

Table  4.5  -  Summary  of  Blowing  Ratio  at  each  Blowing  Level 


Case 

Blowing  Ratio 

Error 

Blowing  0 

0.000128 

±7.18% 

Blowing  1 

0.004640 

±4.18% 

Blowing  2 

0.007695 

±1.72% 

Blowing  3 

0.009853 

±16.9% 

Blowing  4 

0.011552 

±4.32% 

4.2.2  Boundary  Layer  Growth 

Visualization  of  the  boundary  layer  at  each  blowing  level  gave  the  first  indication  that 
transpiration  cooling  had  an  effect  on  boundary  layer  thickness.  Figure  4.5  shows  a  series  of 
shadowgraphs  taken  at  approximately  the  same  time  for  each  blowing  level.  Notice  that  the 
boundary  layer  on  the  non  blowing  wall  is  unchanged  from  the  non  blowing  case  as  it  should  be 
(compare  with  Figure  4.1c),  but  the  increase  in  boundary  layer  thickness  on  the  blowing  side  with 
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(c)  Blowing  2 


(e)  Blowing  4 


Figure  4.5  -  Shadowgraphs  Showing  Boundary  Layer  at  Each  Blowing  Level 
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increasing  blowing  is  clearly  visible.  All  other  visible  features  of  the  flow  including  mach  lines  and 
local  turbulence  patterns  appear  to  be  unchanged  in  each  case  indicating  that  the  disturbance  due 
to  transpiration  was  limited  to  the  boundary  layer. 

As  mentioned,  boundary  layer  growth  was  evaluated  as  a  ratio  of  boundary  layer  thickness 
to  exit  height,  and  this  evaluation  was  accomplished  considering  only  one  side  and  half  the  exit 
height.  The  values  would  be  the  same  if  it  were  assumed  that  both  sides  were  transpiration  cooled 
and  the  entire  exit  height  was  considered.  Table  4.6  gives  the  measured  exit  boundary  layer 
thicknesses,  and  Figure  4.6  shows  the  boundary  layer  growth  data  in  terms  of  comparison  with 
percent  growth. 


Table  4.6  -  Summary  of  Boundary  Layer  Growth  Data 


Case 

Exit  b.l.  thickness, 

Error 

t/(He/2) 

Growth 

Non  Blowing 

1.68  mm  (0.066") 

3.2% 

0.103 

0.0% 

Blowing  1 

1.93  mm  (0.076") 

3.6% 

0.119 

15.5% 

Blowing  2 

2.16  mm  (0.085") 

2.8% 

0.133 

29.1% 

Blowing  3 

2.29  mm  (0.090") 

4.0% 

0.141 

36.8% 

Blowing  4 

2.46  mm  (0.097") 

2.7% 

0.152 

47.5% 

A  quadratic  curve  was  fit  to  this  data  because  this  could  be  useful  for  interpolation 
between  the  data  points  at  this  range  of  blowing  ratios.  It  is  unknown  if  this  curve  fit  would  apply 
to  data  outside  this  range  of  blowing  ratios,  but  is  seems  unlikely  that  the  boundary  layer 
thickness  would  continue  to  increase  at  its  current  rate  for  long  without  leveling  off 

Boundary  layer  thickness  growth  data  along  the  nozzle  did  not  show  the  jump  in  thickness 
at  the  beginning  of  the  porous  material  reported  by  Goldstein  et.al.,  but  there  is  a  good  reason  for 
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this.  Goldstein’s  experiment  involved  a  flat  plate  with  no  pressure  gradient  so  there  was  constant 
blowing  over  the  entire  porous  section.  In  that  case,  when  the  flow  interacted  with  the  injection 
there  was  a  significant  jump  in  boundary  layer  thickness  (See  Figure  1.3).  Because  of  the 
pressure  gradient  present  in  this  experiment,  the  blowing  ratio  starts  small  and  increases  along  the 
porous  area.  This  causes  a  smooth,  blended  growth  along  the  length  of  the  nozzle,  and  this  type 
of  disturbance  in  the  boundary  layer,  as  opposed  to  the  abrupt  disturbance  caused  by  film  cooling 
or  cases  with  no  pressure  gradient,  strengthens  the  case  for  transpiration  cooling. 


- 1 - 1 - ' 

0  0.005  0.01  0.015 

Blowing  Ratio 

Figure  4.6  -  Blowing  Ratio  vs  Boundary  Layer  Growth 


Goldstein’s  lowest  blowing  ratio  (0.0078)  was  the  most  similar  to  those  used  in  this 
experiment,  and  he  showed  approximately  a  60%  increase  in  boundary  layer  thickness,  which  is 
larger  than  the  growth  seen  here  (see  Figure  4.6).  This  is  immediately  seen  as  a  benefit  of 
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transpiration  cooling.  The  smooth  transition  can  be  seen  from  the  shadowgraphs  (See  Figure 
4.5),  and  the  data  in  Figure  4.7. 
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Figure  4.7  -  Boundary  Layer  Growth  vs  Distance  from  Throat 


4.2.3  Wall  Mach  Number 

The  Mach  number  along  the  wall  was  measured  based  on  static  pressure  for  several  runs 
at  each  blowing  level  to  determine  if  blowing  had  a  noticeable  effect  on  wall  Mach  number.  There 
were  small  differences  observed  when  comparing  the  Mach  numbers  on  the  two  sides,  and  several 
factors  such  as  roughness  and  the  laminar  sublayer  thickness  could  have  contributed  to  this 
difference.  Of  more  interest  was  the  fact  that  the  difference  was  the  same  at  every  blowing  level. 
All  wall  Mach  number  data  points  collected  were  plotted  in  Figures  4.8,  but  the  data  was  so 
constant  that  several  data  points  fell  on  top  of  one  another.  The  same  was  true  when  the  Mach 
number  data  was  converted  to  pressure  ratios  in  Figure  4.9.  Therefore,  it  can  be  seen  that 
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blowing  had  no  effect  on  wall  Mach  number  at  any  blowing  ratio.  Because  blowing  did  not  affect 
the  wall  static  pressure,  assuming  a  thin  boundary  layer  for  analysis  of  the  transpired  boundary 
layer  appears  to  remain  a  valid  at  low  blowing  ratios. 


Distance  from  throat  (cm) 

Figure  4.8  -  Mach  Number  vs  Distance  from  Throat 


Also,  the  quality  of  the  model  was  seen  from  the  small  differences  between  the  measured 
static  pressures  and  wall  Mach  numbers  and  the  exact  characteristic  solutions  for  the  nozzle. 
These  differences  were  expected  and  were  attributed  to  the  manufacture  of  the  model.  This  result 
also  shows  that  blowing  at  low  blowing  ratios  does  not  have  major  effects  on  the  flow  field  of  the 
nozzle. 
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Figure  4.9  -  Pressure  Ratio  vs  Distance  from  Throat 


4.3  Performance 
4.3.1  Exit  Mach  Number 

Boundary  layer  growth  with  injection  has  already  been  shown,  which  increased  the  mass 
flow  defect  region  in  the  boundary  layer.  Therefore,  more  of  the  flow  energy  had  to  be  spent  to 
accelerate  the  very  low  velocity  injected  gas  and  overcome  this  defect.  This  loss  of  energy  was 
expected  to  degrade  the  exit  Mach  number  especially  near  the  wall.  Larger  injection  rates  cause 
greater  disturbances  to  the  flow.  When  exit  flow  is  not  uniform  in  a  wind  tunnel,  it  approximates 
freestream  flow  less  accurately,  and  this  introduces  greater  error  into  the  wind  tunnel  data.  For 
transpiration  to  be  an  attractive  method  of  cooling  high  temperature  wind  tunnels,  this  effect  must 
be  quantified. 
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Mach  numbers  were  determined  in  the  same  way  as  the  baseline  case,  but  less  points  per 


case  were  taken.  Data  points  near  the  walls  were  more  important  in  the  blowing  cases  because 
the  effects  were  expected  in  those  regions.  The  centerline  exit  Mach  number  was  also  measured 
as  a  reference  for  each  blowing  level.  Figure  4.10  compares  the  exit  Mach  profile  for  each  case 
with  the  baseline  case  (see  Figure  4.3),  and  Table  4.7  summarizes  the  data. 

The  maximum  error  in  the  Mach  number  data  presented  in  Figure  4.10  was  ±2.87%  so 
there  is  high  confidence  in  the  pitot  probe  measurements.  Errors  in  the  average  non  blowing  and 
blowing  Mach  numbers  were  expected  because  the  Mach  number  did  vary  from  the  centerline  to 
the  wall,  and  with  a  small  number  of  data  points,  the  errors  tended  to  be  large.  This  was  still  a 
better  indication  of  the  effect  of  transpiration  on  each  side  of  the  nozzle  than  a  single  data  point 
near  the  wall. 


Table  4.7  -Summary  of  Exit  Mach  Number  Data 


Case 

Avg 

Error 

Centerline 

Non  Blowing  Avg 

Error 

Blowing  Avg 

Error 

Non  Blowing 

1.886 

±1.06% 

1.92 

1.882 

±1.54% 

1.884 

±2.17% 

Blowing  1 

1.884 

±2.76% 

1.91 

1.900 

_ 

1.870 

±7.43% 

Blowing  2 

1.860 

±1.77% 

1.88 

1.870 

±34.0% 

1.847 

±4.54% 

Blowing  3 

1.862 

±2.42% 

1.90 

1.870 

±34.0% 

1.843 

±6.99% 

Blowing  4 

1.848 

±3.08% 

1.88 

1.875 

±42.4% 

1.820 

±7.63% 
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Figure  4.10a  -  Baseline  and  Blowing  0 


Mach  Number 


Figure  4.10b  -  Baseline  and  Blowing  1 
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Normalized  /distance  from 
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Figure  4.10c  -  Baseline  and  Blowing  2 
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Figure  4.10d  -  Baseline  and  Blowing  3 
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Mach  Number 

Figure  4.10e  -  Baseline  and  Blowing  4 

The  effects  of  transpiration  on  the  exit  Mach  number  were  small  in  this  study,  -2%  in 
overall  average  Mach  number,  and  -3.4%  in  the  blowing  region.  The  average  Mach  number  in  the 
blowing  region  stayed  within  2%  of  the  centerline  Mach  number.  These  results  confirmed  that 
transpiration  would  decrease  the  exit  Mach  number  of  the  nozzle  and  also  showed  that  the  effect 
of  transpiration  was  limited  to  the  blowing  region.  The  decrease  in  exit  Mach  number  had  an 
adverse  effect  on  the  performance  of  the  nozzle  because  Mach  number  is  directly  related  to  exit 
velocity.  In  Figure  4.11,  the  exit  Mach  number  parameters  used  in  this  experiment  are  shown 
relative  to  the  design  Mach  number, 2.0,  against  blowing  ratio  to  visualize  the  effects  of 
transpiration. 
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Figure  4.11  -  Exit  Mach  Number  Trends  vs  Blowing  Ratio 


Blowing  Ratio 

Figure  4.12  -  Decrease  in  Uniformity  w/  Blowing 
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The  other  effect  of  decreased  exit  Mach  number  near  the  walls  (indicated  by  blowing 
average  Mach  number),  was  decreased  uniformity  of  the  exit  flow.  The  importance  of  flow 
uniformity  to  wind  tunnels  has  already  been  stated,  but  it  is  also  a  concern  for  thrust.  Non 
uniform  exit  flow  has  velocity  components  perpendicular  to  the  exit  plane,  which  decrease  thrust. 
Figure  4.12  uses  the  quantity  (1  -  spatial  variance  of  the  exit  flow)  to  show  the  decrease  in 
uniformity  with  increased  blowing. 

4.3.2  Specific  Impulse 

In  rocket  engine  design,  an  important  performance  parameter  is  specific  impulse  or  Isp.  It 
is  a  measure  of  fuel  consumption  efficiency  because  it  is  related  to  mass  flow  rate  and  thrust  by 
Eq.4-1. 

mu,  +AAp,  —p„) 

j  (4-1) 

mgo 

where;  m  =  mass  flow  rate 
Ue  =  exit  velocity 
Ae  =  exit  area 
Pe  =  exit  pressure 
Pa  =  ambient  pressure 

It  was  expected  that  transpiration,  while  adding  mass  flow  to  the  nozzle,  would  degrade 
the  Isp  of  the  nozzle  slightly  because  of  the  exit  Mach  number  decreases  and  the  temperature  of 
the  injected  gas. 
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Specific  impulse  was  calculated  as  an  extension  of  the  blowing  ratio  analysis  with 
BLOWDOT.FOR.  Total  mass  flow  rate  took  into  account  the  mass  flows  due  to  suction  and 
blowing,  In  every  case  except  Blowing  0  the  total  mass  flow  taking  into  account  suction  and 
blowing  was  greater  than  the  freestream  mass  flow  rate.  Calculations  were  originally  made  at  the 
chamber  conditions  specific  to  each  run,  but  the  small  variances  from  run  to  run  caused  the  results 
to  spread  too  much  to  be  useful.  Therefore,  the  data  was  normalized  to  the  average  chamber 
conditions,  and  a  trend  was  noticeable  from  the  results.  Figure  4. 13  shows  the  results. 


Figure  4.13  •  Blowing  Ratio  vs  Isp 


From  Figure  4.13  it  can  be  seen  that  transpiration  had  virtually  no  effect  on  the  specific 
impulse  of  the  nozzle,  and  if  anything,  there  is  a  slight  upward  trend.  The  effect  was  very  small 
due  to  the  transpiration  setup  in  the  nozzle.  First  of  all,  the  blowing  ratios  were  very  small,  which 
means  that  very  httle  mass  was  injected  into  the  flow.  Also,  most  of  the  injection  occurred  near 


63 


the  back  of  the  porous  area  where  the  pressure  drop  was  greatest.  By  the  time  the  freestream  air 
reached  that  point  it  had  cooled  to  a  temperature  that  was  very  close,  and  perhaps  a  little  lower, 
than  the  injected  air.  Therefore,  it  was  not  cooling  the  freestream  air,  which  would  have  been  a 
source  of  performance  loss  in  the  nozzle.  The  combination  of  small  losses  in  exit  velocity  and 
temperature  combined  with  increased  total  mass  flow  and  exit  pressure  acted  to  cancel  each  other 
out  in  this  case.  The  canceling  of  these  trends  resulted  in  the  negligible  effect  on  the  nozzle’s 
performance.  As  an  example  of  the  counteracting  trends.  Figure  4. 14  shows  how  thrust  increased 
slightly  more  rapidly  than  total  mass  flow  rate  at  higher  blowing  ratios. 


Blowing  Ratio 

Figure  4.14  -  Blowing  Ratio  vs  Thrust,  Mass  Flow  Rate,  Exit  Temp,  and  Exit  Pressure 
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V.  Conclusions  and  Recommendations 


The  objectives  of  this  study  were  to  understand  the  differences  between  flow  over  a 
porous  and  non  porous  wall  with  no  blowing,  to  relate  blowing  ratio  to  boundary  layer  thickness 
on  a  porous  wall,  to  relate  blowing  ratio  to  the  uniformity  and  magnitude  of  the  exit  Mach 
number,  and  relate  the  above  results  to  the  performance  of  the  nozzle.  All  of  these  objectives 
were  met,  and  some  additional  conclusions  and  recommendations  were  also  made. 

The  boundary  layer  was  shown  to  be  213%  greater  on  the  blowing  side  of  the  nozzle  than 
on  the  non  blowing  side,  but  this  difference  did  not  affect  the  exit  Mach  number.  The 
conclusion,  then,  is  that  in  this  case,  the  boundary  layer  growth  due  to  the  porous  wall  without 
blowing  was  caused  mostly  by  the  roughness  of  the  material.  Even  without  active  injection, 
however,  suction  and  blowing  along  the  nozzle  due  to  the  pressure  gradient  was  present.  All  of 
these  effects  together  pointed  out  that  the  boundary  layer  over  a  porous  wall  without  active 
injection  was  still  very  complex. 

At  low  blowing  ratios,  the  blowing  side  boundary  layer  thickness  showed  significant 
increases  over  the  non  blowing  case  as  expected.  The  boundary  layer  thickness  increased  by  as 
much  as  47%  over  the  small  range  of  blowing  ratios  studied,  which  made  the  boundary  layer 
thickness  increase  function  steep  at  low  blowing  ratios.  The  increase  along  the  wall  was 
continuous,  however,  which  means  that  the  injection  caused  a  smooth,  blended  flow  disturbance 
along  the  wall,  and  this  is  a  benefit  to  transpiration  cooling  over  film  cooling. 

Exit  Mach  number  was  adversely  affected  by  blowing.  The  average  exit  Mach  number 
decreased  by  2%,  but  more  importantly,  the  average  Mach  number  in  the  blowing  region 
decreased  by  3.4%.  Therefore,  low  blowing  ratios  have  a  noticeable  effect  on  the  exit  flow  of  a 
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nozzle,  and  if  the  entire  nozzle  were  transpiration  cooled,  the  effects  would  be  would  have  been 
greater.  The  Mach  number  profile  across  the  exit  plane  was  also  disrupted  in  terms  of 
uniformity,  as  shown  by  a  50%  increase  in  the  variance  of  the  exit  profile.  Therefore, 
transpiration  cooling  also  has  a  noticeable  effect  on  exit  flow  uniformity  low  blowing  ratios. 

Performance  calculations  showed  that  transpiration  cooling  caused  a  slight  upward  trend 
in  the  specific  impulse  of  the  nozzle.  The  trend,  however,  was  so  slight  that  it  can  be  concluded 
that  at  the  low  blowing  ratios  and  low  thrust  chamber  temperatures  studied  in  this  experiment, 
transpiration  cooling  had  no  effect  on  the  performance  of  the  nozzle  in  terms  of  specific  impulse. 
This  conclusion  is  attributed  to  factors  including  porous  area  geometry  and  the  temperature  range 
of  the  experiment.  If  this  had  been  a  combustion  experiment  with  extremely  high  temperature  in 
the  nozzle  there  almost  certainly  would  have  been  a  decrease  in  performance.  If  an  application 
for  a  low  temperature  nozzle  is  found,  however,  transpiration  cooling  might  prove  useful  for 
more  than  cooling. 

The  geometry  of  the  blowing  area  had  two  effects  on  this  experiment  that  brought  to  light 
one  of  the  major  criticisms  of  transpiration  cooling,  lack  of  flow  control  in  the  presence  of  a 
pressure  gradient.  During  this  experiment,  most  of  the  blowing  occurred  at  the  aft  end  of  the 
porous  area  because  the  flow  tended  to  go  where  the  pressure  was  lowest,  and  this  effect  caused 
three  problems.  First,  blowing  was  concentrated  away  from  the  throat  where  it  was  needed  most. 
Second,  the  exit  boundary  layer  was  probably  thicker  than  it  would  have  been  had  blowing  been 
concentrated  at  the  throat,  and  lastly,  the  temperature  of  the  injected  gas  was  the  same  or  higher 
than  the  freestream  gas,  which  reduced  the  effect  of  injection  on  performance.  Therefore,  unless 
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a  variable  porosity  material  becomes  readily  available,  the  porous  area  in  the  presence  of  a 
pressure  gradient  should  be  such  that  these  effects  are  reduced. 

The  problem  should  not  be  approached  with  wild  abandon,  however,  because  it  is  unclear 
what  combination  of  effects  will  occur  when  the  porous  area  is  concentrated.  The  porous  area 
should  not  be  so  concentrated  as  to  negate  the  transitional  flow  disturbance  benefit  of 
transpiration.  If  this  occurs  then  transpiration  cooling  approximates  film  cooling,  and  this  may 
not  be  desirable.  The  only  certain  thing  about  concentrating  the  porous  area  is  that  that  the 
injection  will  be  applied  at  the  throat  where  it  is  needed.  Therefore,  changing  the  porous  material 
injection  geometry  is  a  logical  extension  of  this  research. 

Two  recommendations  for  future  research  are  made.  First,  more  data  of  the  kind 
presented  here  should  be  taken  at  higher  blowing  ratios  to  extend  the  present  curves  because  it  is 
expected  that  the  trends  established  here  will  continue  at  higher  blowing  ratios.  This  would 
provide  a  more  complete  picture  of  the  effects  of  transpiration  cooling.  Secondly,  the  porous 
area  geometry  should  be  changed  and  a  complete  set  of  curves  be  created  to  compare  with  these. 
This  comparison  should  point  to  the  transpiration  setup  that  provides  the  best  cooling 
effectiveness  with  the  smedlest  performance  loss. 

Heat  transfer  was  shown  in  a  parallel  study  by  Captain  Joseph  Lenertz  to  have  been 
reduced  by  14%  in  the  same  nozzle  used  here,  which  is  significant  considering  the  low  blowing 
ratios  studied  in  this  experiment.  Major  performance  losses  were  not  observed  in  this  experiment 
but  were  proportional  to  the  amount  of  injection,  and  based  on  the  amount  of  heat  transfer 
reduction  observed  in  the  parallel  study,  a  clear  trade  off  exists  between  heat  transfer  reduction 
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and  performance  loss.  Overall,  then,  transpiration  cooling  was  shown  to  have  potential  as  an 
effective  means  of  cooling  with  minimal  performance  loss  in  nozzle  applications. 
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APPENDIX  A 

Gamma  Calculation  for  air/He  Driver 


A.l 


APPENDIX  A  -  Y  Calculation  for  He/Air  Driver  Mixture 


Given:  =  294.43  K 

R,ir  =  287  J/kg-K 
Vine  =  4  kg/mol 
yair=  1.4 
Y//g=  1.667 

Find  Cp  and  c,.  for  air  (Hill  and  Peterson,  p.38): 

Cp!  =  1004.26  J/kg°K 
Cw  =  717.26  J/kg°K 

Find  Cp  and  c,,  for  He: 


8314.3 


J 


R 


kg-K 


He 


kg 

mol 


^p2  ~  ^v2  ^He 

^v2  '  V  He  ~  ^v2  ^He 


C 


v2 


iVne-^) 


2078.5 


J 

kg-K 


2078.5 


J 


Cv2 


*^-^-3116,19^ 


(1.667-1) 
J 


kg-K 


+  2078.5- 


J 


=  5194.69 


kg-K 

Find  densities  of  each  gas  in  the  driver  section: 

p=  pRT 

_ P_ 

^  RT 


J 

kg-K 
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\AApsi  99284,5P^ 

J  1 

(287 - ){29A3K)  (287 - )(294.3^) 

^  ^g  -  ^ 

70.0/757  482500.0Pfl  _  kg 

= - j - - - - - -./88/8  — 

(2078.5^^ — ){29A3K)  (2078.5 - ){29A3K) 

kg -K  kg -K 


Use  density  averaging  method  to  find  mixture  Cp  and  Cv. 

,  P.,A,+PhA.2  (1.17494)(717.26)  +  (.78878)(3116.19)  J 

P...r+PHe  (1.17494+.78878)  '  kg-K 


_  PairCp^  ^Pne^p2  _  (1.1 7494)(1 004.26)  +  (.78878)(5 1 94.69)  J 

Pair^Pn.  (1.1 7494+78878)  ’  kg-K 

/.=  1.598 


=  1006.6 


J 

kg-K 


Ym  calculated  above  was  the  value  used  in  shock  tube  calculations  to  find  chamber 
conditions.  This  calculation  was  accomplished  for  a  range  of  driver  pressure  mixtures  to  create 
the  spreadsheet  shown  on  page  A.  4. 
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Air  Pressure  (psi) 


Spreadsheet  for  calculating  y  over  a  range  of  driver  mixtures 


Air  Properties 
7?  =  287  J/kg-K 
Cp  =  1004.26  J/kg-K 
c,.  =  717,26  J/kg-K 


Helium  Properties 
7?  =  2078.5  J/kg-K 
t>  =  5194.69  J/kg-K 
c,.  =  31 16.19  J/kg-K 


Gas  Temperature  =  294.3  K 


70.00 

70.50 

Helium  Pressure  (psi) 

71.00  71.50 

72.00 

72.50 

14.00 

1.6003 

1.6006 

1,6010 

1.6013 

1.6017 

1.6020 

14.05 

1.6001 

1.6005 

1.6008 

1.6012 

1.6015 

1.6018 

14.10 

1.5999 

1.6003 

1.6006 

1.6010 

1.6013 

1.6017 

14.15 

1.5997 

1.6001 

1.6005 

1.6008 

1.6011 

1.6015 

14.20 

1.5996 

1.5999 

1.6003 

1.6006 

1.6010 

1.6013 

14.25 

1.5994 

1.5997 

1.6001 

1.6005 

1.6008 

1.6011 

14.30 

1.5992 

1.5996 

1.5999 

1.6003 

1.6006 

1.6010 

14.35 

1.5990 

1.5994 

1.5997 

1.6001 

1.6004 

1.6008 

14,40 

1.5989 

1.5992 

1.5996 

1.5999 

1.6003 

1.6006 

14.45 

1.5987 

1,5990 

1.5994 

1.5998 

1.6001 

1.6004 

14.50 

1.5985 

1.5989 

1.5992 

1.5996 

1.5999 

1.6003 

14.55 

1,5983 

1.5987 

1.5990 

1.5994 

1.5998 

1.6001 

14.60 

1.5982 

1.5985 

1.5989 

1.5992 

1.5996 

1.5999 

14.65 

1.5980 

1.5983 

1.5987 

1.5991 

1.5994 

1.5998 

14.70 

1.5978 

1,5982 

1.5985 

1.5989 

1.5992 

1,5996 

14.75 

1.5976 

1.5980 

1.5984 

1,5987 

1.5991 

1,5994 

14.80 

1.5975 

1.5978 

1.5982 

1.5985 

1,5989 

1.5992 

14.85 

1.5973 

1.5976 

1.5980 

1.5984 

1.5987 

1.5991 

14.90 

1.5971 

1.5975 

1.5978 

1.5982 

1.5986 

1.5989 

14.95 

1.5969 

1.5973 

1.5977 

1.5980 

1.5984 

1.5987 

15.00 

1.5968 

1.5971 

1.5975 

1.5979 

1.5982 

1.5986 
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APPENDIX  B 
Computer  Resources 
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Program  SHOCKTUN.FOR 


Purpose:  Use  shock  tube  relations  to  predict  stagnation  conditions  behind  a  shock  wave 

after  reflection  from  a  straight  wall. 

Inputs:  All  inputs  to  SHOCKTUN.FOR  are  from  the  keyboard 

Driver  Length  (ft) 

Driver  Pressure  (psi).  Temperature  (K) 

Driven  Length  (ft) 

Driven  Pressure  (psi).  Temperature  (K) 

Output:  File,  SHOCK. OUT 

SHOCK  TUBE  SUMMARY 
REGION  4  (DRIVER  SECTION) 

Driver  Length(m)  =  1.2I92I 

Gamma  4=  1. 57415 

T4(K)  =  294.430 

P4(atm)  =  4.04082 

REGION  I  (DRIVEN  SECTION) 

Driven  Length(m)  =  5.48647 

Gamma  I  =  1.40000 

TI(K)=  294.430 
PI(atm)=  0.979592 
PRESSURE  RATIO 
p4/pl=  4.12500 
SHOCK  RESULTS 
p2/pl  =  2.29375 

Ms=  1.45222 
Ws(m/s)  =  499.479 

REGION  2  (AFTER  SHOCK) 
up(m/s)  =  u2=  218.866 

T2/T1  =  1.28866 

M2  =  0.560562 

REGION  4  (LEADING  EXPANSION) 

Ue=  -619.077 
a3(m/s)  =  556.246 

REGION  3  (LEADING  EXPANSION) 

Uer=  775.111 

REGION  5  (REFLECTED  SHOCK) 

Mrl=  1.39138 
Wr(m/s)  =  324.386 

p5/p2  =  2.09194 

p5/pl  =  4.79839 

p5/p4=  1.16325 

T5/T2  =  1.24914 

T5/T1  =  1.60971 

STAGNATION  CONDITIONS 
Pt(atm)  =  4.70046 

T5(K)=  473.946 


B.2 


TEST  TIME 
Shock  time  =  10.9844 

Pistion  time  =  25.0677 

Run  time(ms)  =  14.0833 

Time  endvvall  expansion(ms)  =  10.6207 


Program  NOZ2.FOR 

Purpose:  Construct  a  contoured  nozzle  shape  for  a  given  exit  Mach  number  and  geometry  using 
the  method  of  characteristics. 

Inputs:  All  inputs  to  SHOCKTUN.FOR  are  from  the  keyboard 

Desired  number  of  characteristic  reflections  (2  or  4) 
k,  controls  curvature  of  throat  area  (.2  - 1.0  recommeded) 

Exit  height  (input  only  half  the  actual  exit  height  in  inches) 

Ratio  of  specific  heats  (y) 

Desired  exit  Mach  number 

Outputs:  File,  NOZZ  .OUT  -  contains  table  of  characteristic  results,  which  is  that  exact 

solution  at  each  point  in  the  characteristic  lattice. 

Includes:  x,y  coordinates,  Mach  numbers,  Mach  angles 

Char  #  x  y/y,iu.oat  M  0  v 

1  O.OOOOOE+OO  O.IOOOOE+Ol  O.lOOOOE+01  0.34907E-01  O.OOOOOE+00  0.15708E+01 

2  0.29101E+00  0.10051E+01  0.11402E-K)1  0.37771E-01  0.37771E-01  0.10697E+01 

3  0.31491E+00  0.10059E+01  0.11476E+01  0.40635E-01  0.40635E-01  0.10580E+01 

4  0.33881E+00  ().10069E+01  0.11549E+01  0.43499E-01  0.43499E-01  0.10468E+01 

5  0.36272E+00  0.10079E+01  0.11621E-K)1  0.46363E-01  0.46363E-01  0.10362E+01 

6  0.38664E+00  0.10090E+01  0.11692E+01  0.49227E-01  0.49227E-01  0.10261E+01 

7  0.41056E+00  O.lOlOlE+01  0.11762E+01  0.52092E-01  0.52092E-01  0.10163E+01 

8  0,43449E+00  0.10113E+01  0.11831E+01  0.54956E-01  0.54956E-()I  0.1007()E+01 

9  0.45843E+00  0.10126E+01  0.11899E+01  0.57820E-01  0.57820E-01  0.99798E+00 

10  0.48237E+00  0.10140E+01  0.11967E-K)1  0.60684E-01  0.60684E-01  0.98930E+00 

11  0.50632E+00  0.10154E+01  0.12034E+01  0.63548E-01  0.63548E-01  0.98092E+00 

12  0.53028E+00  0.10169E+01  0.12100E+01  0.66412E-01  0.66412E-01  0.97281E+00 

13  0.55425E+00  0.10184E+01  0.12165E+01  0.69277E-01  0.69277E-01  0.96495E+00 

14  0.57823E+00  0.10201E+01  0.12230E-K)1  0.72141E-01  0.72141E-01  0.95733E+00 

15  0.60222E+00  0.10218E+01  0.12294E-K)1  0.75005E-01  0.75005E-01  0.94993E+00 

16  0.62622E+00  0.10235E+01  0.12358E+01  0.77869E-01  0.77869E-01  0.94274E+00 

17  0.65022E+00  0.10254E+01  0.12422E+01  0.80733E-01  0.80733E-01  0.93575E+00 

18  0.67424E+00  0.10273E+01  0.12485E-K)1  0.83598E-01  0.83598E-0I  0.92893E+00 
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File,  WALL.  OUT  -  contains  coordinates  of  upper  wall  and  exact  solution  wall  Mach  numbers 


.X 

O.OOOOE+00 
0, 1087E+00 
0.1176E+00 
().1266E+00 
0.1355E+00 
0.1444E+00 
0.1534E+00 
0.1623E+00 
0.1712E+00 
0.1802E+00 
0.1891E+00 
0,1981E+00 
0.2070E+00 
0.2160E+00 
0.2249E+00 
0.2339E+00 
0.2429E+00 
0.2518E+00 
0.2608E+00 
0.2698E+00 
0.2788E+00 
0.2878E+00 
0.2968E+00 
0.3058E+00 
0.3148E+0() 
0.3238E+00 
0.3328E+00 
0.3418E+00 
0.3508E+00 
0.4499E+00 
0.6864E+00 
0.7140E+00 
0.7413E+00 
0.7683E+00 
0.7950E+00 
0.8214E+00 
0,8477E+00 
0.8739E+00 
0.8999E+00 
0.9258E+00 
0.9517E+00 
0.9776E+00 
0.1003E+01 
0.1029E+01 
0.1055E+01 
0.1081E+01 
0.1107E+01 
0.1133E+01 
0.1159E+01 
0.1185E+01 


y 

0.3735E+00 
0.3754E+00 
0.3758E+00 
0.376  lE+OO 
0.3765E+00 
0.3769E+00 
0.3773E+00 
0.3778E+00 
0.3782E+00 
0.3787E+00 
0.3793E+00 
0.3798E+00 
0.3804E+00 
0.3810E+00 
0.3817E+00 
0.3823E+00 
0.3830E+00 
0.3837E+00 
0.3845E+0() 
0.3852E+00 
0.3860E+00 
0.3868E+00 
0.3877E+00 
0.3885E+00 
0.3894E+00 
0.3904E+00 
0.3913E+00 
0.3923E+00 
0.3933E+00 
0.4048E+00 
0.432  lE+00 
0.4353E+00 
0.4384E+00 
0.4416E+00 
0.4446E+00 
0.4477E+00 
0.4507E+00 
0.4538E+00 
0.4568E+00 
0.4598E+00 
0.4628E+00 
0.4658E+00 
0.4687E+00 
0.4717E+00 
0.4747E+00 
0.4777E+00 
0.4807E+00 
0.4837E+00 
0.4867E+00 
0.4897E+00 


M 

O.lOOOE+01 

0.1140E+01 

0.1148E+01 

0.1155E+01 

0.1162E+01 

0.1169E+01 

0.1176E+01 

0.1183E+01 

0.1190E+01 

0.1197E+01 

0.1203E+01 

0.1210E+01 

0.1217E+01 

0.1223E+01 

0.1229E+01 

0.1236E+01 

0.1242E+01 

0.1248E+0I 

0.1253E+01 

0.1261E+01 

0.1267E+01 

0.1273E+01 

0.1279E+01 

0.1285E+01 

0.1292E+01 

0.1298E+01 

0.1304E+01 

0.1309E+01 

0.1315E+01 

0.1315E+01 

0.1467E+01 

0.1478E+01 

0.1489E+01 

0.1500E+01 

0.1511E+01 

0.1522E+01 

0.1533E+01 

0.1545E+01 

0.1556E+01 

0.1567E+01 

0.1578E+01 

0.1589E+01 

0.1600E+01 

0.1611E+01 

0.1622E+01 

0.1633E+01 

0.1644E+01 

0.1656E+01 

0.1667E+01 

0.1678E+01 
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Program  BLCORR.FOR 


Inputs:  Type  of  gas  (air,He,H2,N2,02),  from  keyboard 

Chamber  Pressure,  from  keyboard 
Chamber  Temperature,  from  keyboard 
File,  WALL, OUT  -  output  from  NOZ2.FOR 

Outputs:  File,  PROPOFX.DAT  -  contains  centerline  data  including  Mach  number,  static 

pressure,  temperature,  density,  flow  velocity,  and  Reynolds  number 


x(in) 

y(in) 

Mach  p(psi)  T(K) 

rho  u(in/s) 

Rex 

Res 

0,000E+00 

0.368E+00 

1.00 

41. 

408,9 

2.42 

405.3 

O.OOE+00 

O.OOE+00 

0.128E+00 

0.370E+00 

1.14 

35. 

389,4 

2.14 

450.9 

0. 14E+06 

0.14E+06 

0.139E+00 

0,370E+00 

1.15 

34. 

388,3 

2.13 

453.4 

0.15E+06 

0.15E+06 

0.150E+00 

0.371E+00 

1.15 

34. 

387.3 

2.12 

455.6 

0.16E+06 

0.16E+06 

0. 160E+00 

0.371E+00 

1.16 

34. 

386.3 

2,10 

457.8 

0.18E+06 

0.18E+06 

0.171E+00 

0.372E+00 

1.17 

34. 

385.3 

2.09 

460.0 

0.19E+06 

0.19E+06 

0.181E+0() 

0.372E+00 

1.18 

33. 

384.3 

2.07 

462.1 

0.20E+06 

0.20E+06 

0.192E+00 

0.373E+00 

1.18 

33. 

383.3 

2.06 

464.3 

0.21E-K)6 

0.21E+06 

0.202E+00 

0.373E+00 

1.19 

33. 

382.4 

2.05 

466.4 

0.22E+06 

0.22E+06 

0.213E+00 

0.374E+00 

1.20 

32. 

381.4 

2.03 

468.6 

0.23E+06 

0,23E+06 

0.223E+00 

0,375E+00 

1.20 

32. 

380.5 

2.02 

470.4 

0.25E+06 

0.25E+06 

0.234E+00 

0,375E+00 

1.21 

32. 

379.5 

2,01 

472.5 

0.26E+06 

0.26E+06 

0,245E+()0 

0,376E+00 

1.22 

32. 

378.5 

2.00 

474.6 

0.27E+06 

0.27E+06 

0.255E+00 

0.377E+00 

1.22 

31. 

377.7 

1.99 

476.4 

0.28E+06 

0.28E+06 

0.266E+00 

0.377E+00 

1.23 

31. 

376.8 

1.97 

478.2 

0.29E+06 

0.29E+06 

0.276E+00 

0.378E+0() 

1.24 

31. 

375.8 

1.96 

480.3 

0.30E+06 

0.30E+06 

0.287E+00 

0,379E+0() 

1.24 

30. 

375.0 

1.95 

482.1 

0.31E+06 

0.31E+06 

0.298E+00 

0.380E+00 

1.25 

30. 

374.1 

1.94 

483.8 

0.33E+06 

0.33E+06 

0.308E+00 

0.381E+00 

1,25 

30. 

373.1 

1.93 

485.9 

0.34E+06 

0.34E+06 

0.319E+00 

0.382E+00 

1.26 

30. 

372,3 

1,92 

487.7 

0.35E+06 

0.35E+06 

File,  WALLC.OUT  -  contains  corrected  wall  coordinates 


x(in) 

O.OOOOE+00 
0.1285E+00 
0.1390E+00 
0.1496E+00 
0.1601E+00 
0.1707E+00 
0.1813E+00 
0.1918E+00 
0.2024E+00 
0.2130E+00 
0.2235E+00 
0.234  lE+00 
0.2447E+00 
0.2553E+00 
0.2659E+00 
0.2765E+00 


yorg(in) 

0.3679E+00 

0.3701E+00 

0.3705E+00 

0.3709E+00 

0.3714E+00 

0.3719E+00 

0.3724E+00 

0.3729E+00 

0,3735E+00 

0.3741E+00 

0,3747E+00 

0.3753E+00 

0.3760E+00 

0.3768E+00 

0,3775E+00 

0.3783E+00 


ycor(in) 

0.3'679E+00 

0.3712E+00 

0.3717E+00 

0.3722E+00 

0.3727E+00 

0.3733E+00 

0,3739E+00 

0.3745E+00 

0.3751E+00 

0,3758E+00 

0.3765E+00 

0.3771E+00 

0.3779E+00 

0.3787E+00 

0.3795E+00 

0.3804E+00 


d*x(in)  d*s(in) 
O.OOOOE+00  O.OOOOE+00 
0.5628E-03  0.5628E-03 
0.5992E-03  0.5993E-03 
0.6355E-03  0.6356E-03 
0.6709E-03  0.6710E-03 
0.7062E-03  0.7064E-03 
0.7410E-03  0.7413E-03 
0.7752E-03  0.7754E-03 
0.8092E-03  0.8096E-03 
0.8430E-03  0.8434E-03 
0,876IE-03  0.8763E-03 
0.9092E-03  0.9097E-03 
0.9420E-03  0.9426E-03 
0.9746E-03  0,9752E-03 
0.1007E-02  0.1008E-02 
0.1039E-02  0.1040E-02 


B.5 


0.2871E+00  0.3791E+00  0.3812E+00  0.1071E-02  0.1072E-02 
0.2977E+00  0.3799E+00  0.3821E+00  0,1102E-()2  0.1103E-02 
0.3083E+00  0.3808E+00  0.383 lE+00  0.1134E-02  0.1135E-02 
0.3189E+00  0.3817E+00  0.3840E+00  0.1165E-02  l).1166E-02 
0.3295E+00  ().3827E+00  0.385 lE+OO  0.1196E-02  0.1197E-02 
0.3401E+00  0.3836E+00  0.3861E+00  0.1227E-02  0.1228E-02 
0.3508E+00  0.3846E+00  0.3871E+00  0.1258E-02  0.1259E-02 
0.36I4E+00  0.3857E+00  0.3883E+00  0.1288E-02  0.1290E-02 
0.3720E+00  0.3867E+00  0.3893E+00  0.1318E-02  0.1320E-02 
0.3827E+00  0.3878E+00  0.3905E+00  0.1349E-02  0.1351E-02 
0.3933E+00  0.3889E+00  0.3917E+00  0.1379E-02  0.1381E-02 
0.4040E+00  0.390 lE+00  0.3929E+00  0.1409E-02  0.1411E-02 
0.4147E+00  0.3913E+00  0.3942E+00  0.1439E-02  0.1441E-02 
0.4341E+00  0.3935E+00  0.3965E+00  0.1492E-02  0.1495E-02 
0.6960E+00  0.4238E+00  0.4282E+00  0.2190E-02  0.2197E-02 
0.7256E+00  0.4272E+00  0.4317E+00  0.2265E-02  0.2273E-02 
0.7548E+00  0.4306E+00  0.4353E+00  0.2339E-02  0.2347E-02 


Program  WALL.FOR 

Inputs:  File,  WALLC  OUT  -  output  from  BLCORR.FOR 

Outputs:  File,  NOZ.TOP  -  contains  x,y  coordinates  of  subsonic  and  supersonic  sections 


X  y 


O.OOOOE+00 

0.4184E+01  -  First  line  of  output  data  is  beginning  of  converging  section 

0.2000E+01 

0.4184E+01 

0.2092E+01 

0.4179E+01 

0,2184E+01 

0.4165E+01 

0.2276E+01 

0.4143E+01 

0.2367E+01 

0.4112E+01 

0.2459E+01 

0.4073E+01 

0.2551E+01 

0.4026E+01 

0.2643E+01 

0.3973E+01 

0.2735E+01 

0.3912E+01 

0.2827E+01 

0.3845E+01 

0.2918E+01 

0.3772E+01 

0.3010E+01 

0.3693E+01 

0.3102E+01 

0.3609E+01 

0.3194E+01 

0.3521E+01 

0.3286E+01 

0.3427E+01 

0.3378E-K)1 

0.3330E+01 

0.3469E+01 

0.3229E+01 

0.3561E+01 

0.3125E+01 

0.3653E+01 

0.3017E+01 

0.3745E+01 

0.2908E+01 

0.3837E+01 

0.2796E+01 

0.3929E+01 

0.2682E+01 

0.4020E+01 

0.2567E+01 

0.4112E+01 

0.245  lE+01 

0.4204E+01 

0.2334E+01 
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0.4296E+01 

0.2218E+01 

0.4388E+01 

0.2101E+01 

().4480E+01 

0.1985E+01 

0.457  lE+01 

0.1870E+01 

0.4663E+01 

0.1756E+01 

0.4755E+01 

0.1644E+01 

0.4847E+01 

0.1534E+01 

0.4939E+01 

0.1427E+01 

0.503  lE+01 

0.1323E+01 

0.5122E+01 

0.1222E+01 

0.5214E+01 

0.1124E+01 

0.5306E+01 

0.1031E+01 

0.5398E+01 

0.9424E+00 

0.5490E+01 

0.8585E+00 

0.5582E+01 

0.7798E+00 

0.5673E+01 

0.7068E+00 

0.5765E+01 

0.6398E+00 

0.5857E+01 

0.5793E+00 

0.5949E+01 

0.5255E+00 

0.604  lE+01 

0.4790E+00 

0.6133E+01 

0.4400E+00 

0.6224E+01 

0.409  lE+00 

0.6316E+01 

0.3865E+00 

0.6408E+01 

0.3726E+00 

0.6500E+01 

0.3679E+00  -  Throat  -  Previously  corrected  diverging  section  follows  from  here 

0.6639E+01 

0.3717E+00 

0.6660E+01 

0.3727E+00 

0.6681E+01 

0.3739E+00 

0.6702E+01 

0.375  lE+00 

0.6723E+01 

0.3765E+00 

0.6745E+01 

0.3779E+00 

0.6766E+01 

0.3795E+00 

0.6787E+01 

0.3812E+00 

0.6808E+01 

0.383  lE+00 

File,  NOZ.BOT  -  contains  x,y  coordinates  of  subsonic  and  supersonic  sections  if  bottom 
wall  was  straight 

Output  has  same  general  format  with  slight  differences.  File  not  really  usefiil 
unless  constructing  only  one  half  a  nozzle. 


■X 

O.OOOOE+00 

0.2000E+01 

0.2092E+01 

0.2184E+01 

0.2276E+01 

0.2367E+01 

0.2459E+01 

0.2551E+01 

0.2643E+01 

0.2735E+01 


y 

-0.3816E+01 
-0.3816E+01 
-0.381  lE+01 
-0.3797E+01 
-0.3775E+01 
-0.3744E+01 
-0.3705E+01 
-0.3658E+01 
-0.3605E+01 
-0.3544E+01 
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Program  BLOWDOT.FOR 


Inputs:  File,  BLOWIN.DAT 

HE8  1.39  70.71  485.0918  14.276  71.7  29.067  1.886 
HE9  1.39  70.71  485.0918  14.276  71.8  29.067  1.886 
HEIO  1.39  69.255  475.1057  14.276  71.6  29.067  1.886 
TlMEl  1.39  69.983  480.0988  14.295  71.7  29.105  1.886 
T1ME2  1.39  69.983  480.0988  14.295  71.7  29.105  1.886 
T1ME3  1.39  69.983  480.0988  14.295  71.7  29.105  1.886 
BLOW2  1.39  69.983  480.0988  37.352  290.15  29.2  1.884 
BLOW4  1.39  69.983  480.0988  37.617  293.261 1  29.2  1.884 
BLOW5  1.39  69.983  480.0988  37.737  290.65  28.985  1.884 
BLOW6  1.39  69.983  480.0988  37.799  292.0944  28.988  1.884 
1M1D2  1.39  69.983  480.0988  45.431  288.7056  29.125  1.86 
1MID3  1.39  69.983  480.0988  45.513  288.7056  29.125  1.86 
1M1D4  1.39  69.983  480.0988  45.529  289.261 1  29. 12  1.86 
1MID5  1.39  69.983  480.0988  45.658  288. 15  29. 12  1.86 
1MID6  1.39  69.983  480.0988  45.721  288.2611  29.12  1.86 
2M1D1  1.39  69.983  480.0988  50.805  287.5944  29.12  1.862 
2M1D2  1.39  69.983  480.0988  50.869  287.8722  29.12  1.862 
MAXI  1.39  69.983  480.0988  55.44  289.8722  29.037  1.848 
MAX2  1.39  69.983  480.0988  55.352  287.7056  29.045  1.848 
MAX3  1.39  69.983  480.0988  55.042  288.0944  29.045  1.848 
MAX4  1.39  69.983  480.0988  55.607  288.2611  29.045  1.848 
MAX5  1.39  69.983  480.0988  49.996  287.3722  29.045  1.848 

Format:File  Name,  y,  chamber  pressure  (psi),  chamber  temperature  (K),  plenum  pressure  (psi), 
plenum  temperature  (F  of  K),  atmospheric  pressure  (in  Hg),  exit  Mach  number 

Outputs:  File,  BLOWRAT.OUT 


Test  Conditions: 

Run:HE8 

Gamma  =  1.3900 
Pa  =  14.2764  psi 
PO  =70.7100  psi 
TO  =485.0918  K 
Ppl  =  14.2760  psi 
Tpl  =295.2056  K 

Initial  Calculations: 
mdotinf=  1.6946  kg/sec 
RHOUinf  =892.4648  kg/m^-s 
rhoinj  =  1.1618  kg/m’ 


Heat  Flux  Gauge  Blowing  Ratios: 

B(l)  =0.000000 
B(2)  =0.000000 
B(3)  =  0.000028 
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Total  Blowing  Ratio  Calculations: 
Isum  =  1.099075  m 
RHOUsum  =  0.000144  kg/m' -s 
BT  =0.000131 

Performance  Calculations: 

Pos  mdot  =  0.000157  kg/sec 
Neg  mdot  =  -.003254  kg/sec 
Tot  mdot  1.6915  kg/sec 
Pe  =  10.7833  psi 

Te  =  286.2008  K 

lie  =  637.2729  m/s 

Thnist=  998.48  N 
Isp  =60.1711  sec 


))(  s((  9|c  9(c 
9)c  :t(  :1c)|c 
)|c  :)c  :|c 


*****  Program  BLOWDOT  calculates  the  total  mass  flow  through  the  porous 
*****  material  for  given  test  conditions.  Blowing  ratios  are  calculated  in  one  manner 
*****  which  will  be  indicated  for  heat  flux  purposes,  and  another  way  (the  Keener  #, 

*****  which  will  be  indicated)  for  performance  purposes.  The  program  performs  a 

j((  3|<)|( 

*****  summation  of  surface  length  and  mass  flow  rate  along  the  porous  surface  at  intervals***** 

*****  defined  by  the  coordinates  of  the  model  in  the  blowing  ratio  calculation.  *  *  *  *  * 

*****  ***** 

*****  Author:  2Lt  David  N.  Keener  513-255-3636  xl249  ***** 

*****  September  1994  ***** 

****************************************************************************** 

*** 

PROGRAM  BLOWDOT 


IMPLICIT  NONE 


REAL*8  rhoinj,Ppl,R,Tpl,uinj,ainj,SA,mdot,mdotinf,B(3),BT,Knum 
REAL*8  mdotsum,CF,P,Astar,in2msq,RHOUinf,PO,TO,psi2Pa,linj 
REAL*8  M(54),x(54),y(54),GAMMA,xsa,Aexit,RHOUinj, RHOUsum, Isum 


INTEGERM  i 


CHARACTER*  8  NAME 

OPEN(  1 , 'c:\qed\data\isentrop.daf) 
OPEN(2,'c:\qed\data\blowin.dat') 
OPEN(3 ,  'c :  \qed\data\blo  wrat .  ouf ) 
OPEN(4,'c  :\qed\data\plotmdot .  ouf) 
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*****  Variable  and  Output  File  Initialization  ***** 

write(3,*)  'Blowing  Ratios  at  each  heat  flux  gauge  and  Total  Bio 
owing  Ratio' 
write(3,*) 

write(3,*)  'Test  Conditions:' 
write(4,*)  'Mass  flow  data  at  each  interval' 
write(4,*) 

write(4,*)  'Test  Conditions:' 

SA=0.0D0 
nidot=0.0D0 
mdotsum=0.0D0 
lsum=0.0D0 
RHOUsuni=O.ODO 
DO  1=1,3 
B(i)  =  O.ODO 
ENDDO 
BT=0.0D0 

*****  Constants  ***** 

R  =  287.0D0 
in2msq  =  .000645 16D0 
psi2Pa  =  6894,757200 
Astar  =  .367900*8, 0D0*in2msq 
Aexit  =  .6395D0*8.0D0*in2msq 

*****  Read  in  data  needed  for  problem  ***** 

DO  i=l,54 

READ(l,*)x(i),y(i),M(i) 

ENDDO 

READ(2,*)  NAME,GAMMA,PO,TO,Ppl,Tpl  Ipressure  in  psi! 

*****  Convert  plenum  temperature  to  Kelvin  if  it  is  not  input  that  way 
IF  (Tpl  It.  lOO.ODO)  THEN 
Tpl  =  (Tpl+459.67D0)/1.8D0 
ENDIF 

write(3,*)  'Run:',NAME 
write(3,*) 

write(3,30)  'Gamma  =', GAMMA 
write(3,31)'P0  =',P0 

write(3,32)  'TO  =',T0 
write(3,31)'Ppl  =',Ppl 
write(3,32)  'Tpl  -,Tpl 
write(4,30)  'Gamma  GAMMA 
write(4,31)'P0  =',P0 


!J/(kg-K)! 

!m^2! 

!m^2! 


write(4,32)  'TO  =',TO 
write(4,3 1 )  'Ppl  =',Ppl 
write(4,32)  'Tpl  =',Tpl 
PO  =  P0*psi2Pa 
Ppl  =  Ppl*psi2Pa 

*****  Initial  calculations  for  iteration  ***** 
write(3,*) 

write(3,*)  'Initial  Calculations:' 
write(4,*) 

write(4,*)  'Initial  Calculations:' 

nidotinf^DSQRT(GAMMA/(R*T0)*((2.0D0/(GAMMA+1.0D0))** 
c  ((GAMMA+1.0D0)/(GAMMA-1.0D0))))*Astar*P0  !kg/sec! 

write(3,30)  'mdotinf  =',mdotinf 
write(4,30)  'mdotinf  =',  mdotinf 

RHOUinf^mdotinf'Astar  !kg/(m^2-sec)! 

write(3,30)  'RHOUinf  =',RHOUinf 
write(4,30)  'RHOUinf  =', RHOUinf 

rhoinj=Ppl/(R*Tpl)  !kg/m^3! 

write(3,30)  'rhoinj  =',rhoinj 

write(4,30)  'rhoinj  =', rhoinj 

write(3,*) 

write(4,*) 

write(4,40) 

*****  Begin  summation  iteration.  For  this  experiment  only  53  points  were  needed 
*****  along  the  nozzle.  This  can  be  changed  appropriately. 

DO  i=2,53 

*****  Correction  factors  were  fit  to  correct  isentropic  predictions  from  original  design  to 
*****  to  measured  data  to  increase  accuracy  of  blowing  ratios. 

IF  (M(i)  It.  1.69I9D0)  THEN 

CF=0.0155*M(i)**3.0D0  +  0.9722D0 
ENDIF 

IF  (M(i)  .ge.  I  .6919D0)  THEN 

CF=0.0I32*M(i)**6.0D0  -  0.0356*M(i)**4.0DO  +  1.0302 
ENDIF 

P=P0/((I.0D0+(GAMMA-1.0D0)/2.0D0*M(i)**2.0D0)**(GAMMA/ 
c  (GAMMA-1. 0D0)))*CF  !Pa! 

*****  Delta  P  across  porous  material  because  there  are  four  layers  of 
*****  material  between  plenum  and  nozzle  wall 

uinj=(DABS(((Ppl/psi2Pa)-(P/psi2Pa))/4.0D0)/.03243D0)** 
c  (1.0D0/.9585D0)*.3048D0/60.0D0  !m/s! 

IF  ((Ppl  -  P)  It.  O.ODO)  THEN 


Hie  9|e  >(c 
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uinj  =  uinj*(-1.0D0) 

ENDIF 

linj=DSQRT((((x(i+l)+x(i))/2-(x(i)+x(i-l))/2)**2.0D0)+ 
(((y(i+l)+y(i))/2-(y(i)+y(i-l))/2)**2.0D0)) 
ainj=linj*3,0D0*in2msq  !m^2! 

RHOUinj  =  rhoinj*uinj 

mdot=RHOUinj*ainj  !kg/sec! 

IF  (i  eq,  11)  THEN 

xsa  =  (y(i)  +.7976D0  *  (y(i+l)-y(i)))  *  8.0D0  *  in2msq 
B(l)  =  ((mdotsum  +  .2853D0*nidot)/(SA  +  .2853D0*ainj))/ 
(mdotinf/xsa) 

IF  (B(l).lt.  O.ODO)  THEN 
B(1)  =  0.0D0 
ENDIF 

write(*,*)  'Writing  Blowing  Ratio  1  to  file' 

WRITE(3,20)  'B(1)=',B(1) 

ENDIF 

IF  (i  .eq.  23)  THEN 

xsa  =  (y(i+l)  +.7976D0  *(y(i+2)-y(i+I)))  *  8.0D0  * 
in2msq 

B(2)  =  ((mdotsum  +  .7828D0*mdot)/(SA  +  .7828D0*ainj))/ 
(mdotinf/xsa) 

IF  (B(2)  It.  O.ODO)  THEN 
B(2)  =  O.ODO 
ENDIF 

write(*,*)  'Writing  Blowing  Ratio  2  to  file' 

WRITE(3,20)  'B(2)  =',B(2) 

ENDIF 

IF  (i  eq.  47)  THEN 

xsa  =  (y(i)  +.7976D0  *(y(i+l)-y(i)))  *  8.0D0  *  in2msq 
B(3)  =  ((mdotsum  +  .1201D0*mdot)/(SA  +  .1201D0*ainj))/ 
(mdotinf/xsa) 

IF  (B(3)  It.  O.ODO)  THEN 
B(3)  =  O.ODO 
ENDIF 

write(*,*)  'Writing  Blowing  Ratio  3  to  file' 

WRITE(3,20)  'B(3)  =',B(3) 

ENDIF 

lsum=lsum+linj 

RHOU  sum=RHOU  sum+(RHOUinj/ (mdotinf/(y(i)  *  8 .  ODO  *  in2msq)))  * 
linj 

SA=ainj+SA  !m^2! 

mdotsum=mdot+mdotsum  Ikg/sec! 

IF  (mdotsum  .le.  O.ODO)  THEN 
mdotsum  =  O.ODO 


lsum=O.ODO 

RHOUsum=O.ODO 

ENDIF 

write(4, 1 0)  (i- 1  ),x(i),CF,(P/psi2Pa),mdot,mdotsum 
ENDDO 

BT  =  (mdotsum/SA)/(mdotinf/Aexit) 
c  BTavg  =  (mdotsum/SA)/(RHOUsum/lsum) 

Knum  =  RHOUsum/lsum 
WRITE(3,20)  'BT  =',BT 
WRITE(3,50)  'Keener  #  =',Knum 
c  WRITE(3,50)  'Isum  =',lsum 

c  WRITE(3,50)  'RHOUsum  =',RHOUsum 

STOP 

10  format(I2,3X,F6.4,3X,F6.4,3X,F5.2,3X,F10.8,3X,F10.8) 

20  format(A7,F8.6) 

30  format(A9,F8.4) 

31  format(A9,F8.4,' psi') 

32  format(A9,F8.4,'  K') 

40  format(8X,'x',7X,'CF',7X,'P',9X,'mdot',7X,'mdotsum') 

50  format(A10,F10.6) 

END 
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APPENDIX  C 
Calibration  Data 


C.l 


Number 

Serial  Number 

Location 

Slone 

1 

WB80 

Tube  1 

0.15838 

2 

44AM 

Tube  2 

0.14378 

3 

SNJ82P 

Plenum 

0.01836 

4 

TM73 

Non  Blowing  1 

0.06754 

5 

29BA 

Converging  Section 

0.02640 

6 

TM83 

Blowing  1 

0.06895 

7 

TN05 

Blowing  2 

0.05519 

8 

TN04 

Blowing  3 

0.07007 

9 

WL43 

Non  Blowing  2 

0.05073 

10 

39BP 

Non  Blowing  3 

0.05344 

Pitot  0.11203 


Very  high  confidence  was  placed  on  these  calibrations  because  the  maximum  error  from  a 
linear  best  fit  in  any  of  these  transducers  was  0.485%  of  the  standard  deviation  in  the  data. 
Therefore,  it  is  unlikely  that  error  in  any  of  the  pressure  data  presented  here  was  due  to  calibration 
error. 
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APPENDIX  D 
Procedures 


D.l 


Shock  Tube  Run  Procedure  -  Helium  Driver  w/  Shadowgraph 

**  NOTE:  If  running  open-air  shock,  make  sure  EVERYONE  IN  ROOM  IS 
WEARING  EAR  PROTECTION!!!! 

**  NOTE:  When  taking  shadowgraphs,  as  many  possible  light  sources  must  be  eliminated 
from  the  room  including  computer  monitors. 

**  NOTE:  Before  starting  a  test  including  shadowgraphs,  turn  out  all  room  lights  except 
those  controlled  by  the  switch  near  the  tube  controls. 

1 .  With  tube  closed,  open  air  valve  and  “P4.” 

2.  Allow  house  air  to  run  through  the  tube  for  5  minutes  to  purge  any  excess  helium  from 
previous  runs  or  any  stagnant  ambient  air. 

During  the  purge  period: 

3.  Load  film  cartridge  into  film  holder. 

4.  Load  film  holder  into  camera. 

5.  Reset  Model  436  Proportional  Delay  Generator  and  Model  453  Delay  Generator. 
When  purge  period  is  complete: 

6.  Close  air  valve  and  “P4.” 

7.  Open  tube  and  remove  spent  mylar  diaphragm. 

8.  Reset  plunger. 

9.  Close  and  secure  tube  with  hydraulic  jack. 

10.  Open  helium  bottle  and  ball  valve. 

1 1 .  Select  “One  Shot”  on  Nicolet  system. 

12.  Turn  down  brightness  on  computer  monitor. 

13.  Open  “Plunger”  valve  and  verify  pressure  to  plunger. 

14.  Open  “P4  Gauge.” 

15.  Slowly  open  “P4.”  P4  must  be  increasing  at  a  slow  rate. 
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16.  Turn  out  remaining  room  lights. 

17.  Expose  film. 

18.  Allow  pressure  to  rise  in  driver  to  desired  pressure  then  close  “P4.” 

19.  Repeatedly  open  and  close  “P4”  until  driver  settles  at  desired  pressure. 

20.  Verify  no  personnel  standing  near  exit  of  shock  tube  and  PhD  students  in  room  are 
aware  of 

impending  test. 

21.  Initiate  test  by  depressing  the  plunger  arm. 

22.  Relieve  excess  plunger  pressure  by  opening  “Plunger  Bleed.” 

23.  Relieve  excess  driver  pressure  by  opening  “P4  Bleed.” 

24.  Close  helium  bottle  and  ball  valve. 

25.  Close  film  cartridge  and  begin  image  developing. 

26.  Turn  on  room  lights. 
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Shock  Tube  Procedure  -  Air  Only  Driver 

**  NOTE:  If  mnning  open-air  shock,  make  sure  EVERYONE  IN  ROOM  IS 
WEARING  EAR  PROTECTION!!!! 

1.  Close  “Plunger”  valve  and  “P4”  valve. 

2.  Open  “Plunger  Bleed”  and  “P4  Bleed,” 

3.  Reset  plunger. 

4.  Close  “Plunger  Bleed”  and  “P4  Bleed.” 

5.  Load  shock  tube  with  appropriate  thickness  mylar  diaphragm. 

6.  Open  “Plunger”  valve  and  verify  pressure  to  plunger. 

7.  Open  “P4  Gauge.” 

8.  Slowly  open  “P4.” 

9.  Select  “One  Shot”  on  Nicolet  system. 

10.  Allow  pressure  to  rise  in  driver  to  desired  pressure  then  close  “P4.” 

1 1 .  Open  and  close  “P4”  until  driver  settles  at  desired  pressure. 

12.  Verify  no  personnel  standing  near  exit  of  shock  tube  and  PhD  students  in  room  are 
aware  of  impending  test. 

13.  Run  shock  by  depressing  plunger  arm. 

14.  Open  shock  tube  and  remove  spent  mylar. 
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Throat  Boundary  Layer  Thickness  Calculation 


E.l 


Discharge  Coefficient: 


Q=l- 


'  -2.128 


r^+iY  _ 

I  2  J  I  /+1 


-  + 3.266  r/?  ^+0.9428 


(r-i)(x  +  2) 

ir+\y^ 


(2-13) 


4  R 

R=R.\- 


(2-14) 


where:  /=  ratio  of  specific  heats  (Cp/C\)  =1.39 
Re*  =  Reynolds  number  at  the  throat 
R*  =  throat  radius  =  9.34466x10''’  m 
Rc  =  radius  of  curvature  of  the  throat  =  0.9"  =  .02286  m 
R  =  constant  used  to  calculate  discharge  coefficient 

The  only  other  quantity  necessary  to  solve  Eq.2-14  is  Re*,  which  is  given  by  Eq.E-1. 


Re*  = 


♦  *  7-^• 

p  u  D 
* 


(E-1) 


where:  p*  =  density  of  gas  at  the  throat  (kg/m’) 
u*  =  velocity  of  gas  at  the  throat  (m/s) 

D*  =  diameter  of  the  throat  (m)  =  1.868932x10'^  m 

p*  =  viscosity  at  the  throat  (N-s/m^)  =  229.4x10'^  N-s/m^  (Hill  and  Peterson,  694) 


Representative  chamber  conditions  were  taken  to  be:  477  kPa  (69.24  psi) 

475  K 

By  isentropic  relations,  throat  conditions  were:  253  kPa  (36.69  psi) 

397.49  K 
E.2 


Therefore,  because 


P 

* 

u 


P 

RT* 


h  K  T') 


at  the  throat, 

p  =2.21718  kg/m3 
u*  -  398.208  m/s 
and 


Re*  =  719,303.615. 

It  follows,  then,  from  Eqs.  2-13  and  2-14  that  R  =  4.83605x10^  and  Co  =  .996097.  Using 
Eq.2-15,  6t  is  calculated  to  be  0.0182  mm  (0.0007179"). 


1-Cp 

R*  2 


(2-15) 


This  value  was  calculated  because  the  boundary  layer  was  not  measurable  near  the  throat 
from  shadowgraphs  because  it  was  so  thin  there.  This  value  was  used  as  a  data  point  in  this 
report  when  boundary  layer  at  the  throat  was  plotted  or  quoted  because  it  was  known  that  the 
boundary  layer  did  exist  at  the  throat,  and  this  is  a  good  estimate  of  its  thickness. 
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APPENDIX  F 


Selecting  a  Porous  Material 


F.l 


The  porous  material  used  in  this  experiment  was  chosen  based  on  strength  mass  flow  rate 
considerations.  The  pressure  the  material  was  required  to  hold  was  the  major  factor  in  choosing 
the  type  of  metal,  used,  and  the  desired  mass  flow  rate  through  the  material  determined  the 
porosity  chosen.  Mott  Metallurgical  Company  was  identified  as  the  source  of  the  material,  and 
analysis  was  accomplished  based  on  strength  and  mass  flow  data  provided  by  Mott. 

Strength  Analysis 

Design  plenum  pressure:  552  kPa  (80  psi) 

Maximum  plenum  pressure:  690  kPa  (100  psi)  (defined  by  the  limit  of  the  pressure  transducers) 
Blowing  area:  58  cm^ 

Maximum  pressure  force  on  the  porous  material:  4  kN 
Maximum  pressure  on  porous  area:  690  kPa 

Pressure  limit  of  3 16L  sintered  stainless  steel:  53MPa 

The  316L  sintered  stainless  steel  was  the  least  strong  material  available,  and  the  pressure 
shown  above  was  material  with  the  greatest  porosity  available.  Therefore,  316L  stainless  steel 
was  more  than  adequate  to  withstand  the  expected  pressure  force,  and  this  material  was  selected 
as  the  type  of  porous  material,  but  porosity  grade  also  had  to  be  selected. 
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Mass  Flow  Analysis 

The  porosity  grade  depended  on  the  desired  mass  flow  through  the  porous  material.  Mott 
Metallurgical  manufactured  their  material  in  1.57  mm  (.062")  thick  sheets,  and  in  order  to  mount 
pressure  transducers  and  heat  flux  gauges  properly  four  sheets  of  material  were  used.  Mass  flow 
through  porous  material  is  related  to  the  pressure  drop  across  the  material,  and  since  all  mass  flow 
data  provided  by  Mott  was  quoted  for  one  layer,  all  pressure  drops  in  the  following  calculations 
had  to  be  divided  by  four  to  take  each  layer  into  account.  Arbitrarily,  the  desired  blowing  ratios 
were  chosen  to  be  0.01  and  0.05  at  the  throat. 

Predicted  chamber  conditions:  538  kPa  (78  psi) 

491  K 

and  isentropic  relations  allowed  calculation  of  pu  (where  *  indicates  the  throat). 

The  injection  parameter,  is  then  calculated  by  Eq.F-1. 

p^u^  =  B-{pV)  (F-1) 

where:  B  =  blowing  ratio 

p  =  gas  density  (kg/m^) 
u  =  gas  velocity  (m/s) 

Pressure  to  the  plenum  was  provided  with  bottled  air. 

Assumed  bottle  conditions:  13.8  MPa  (2000  psi) 

297  K 

Plenum  temperature  (from  Mollier  diagram  for  nitrogen  at  constant  enthalpy):  272  K 
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Using  the  estimated  plenum  temperature,  the  plenum  pressure,  and  the  perfect  gas  law,  pp 
was  calculated.  Assuming  that  p,  =pp  because  the  porous  material  was  thin,  was  calculated. 

From  a  log.-log  chart  relating  mass  flow  to  pressure  drop  provided  by  Mott  (see  Fig  F.  1), 
equations  of  the  form  (F-2)  were  created. 

(F-2) 

where;  Ap  =  pressure  drop  (psi) 

A,B  =  constants  determined  from  the  chart 
Ui  =  injection  velocity  (ft/min) 

Inserting  the  design  pressure  drop  and  injection  velocity  in  the  correct  units  into  Eq.F-2 
determined  that  2  micron  porosity  was  suited  for  B=0.01  and  10  micron  porosity  was  suited  for 
B=0.05.  In  reality,  only  the  2  micron  porosity  material  was  used,  and  this  analysis  was  validated 
to  within  2%  of  Mott  data. 


I  e  lo  '  so  100  BOO  1000 

AIR  FLOW,  CFM/FT’ 


Figure  F.l  -  Mass  Flow  vs  Pressure  Drop  (Mott) 
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